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Preface 


The  purpose  of  this  study  was  to  design  a  digital  communication 
model  and  specific  jamming  techniques  which,  together  or  separately, 
would  be  used  to  test  the  performance  of  various  error  correcting 
codes.  Much  of  my  research  was  spent  becoming  familiar  with  current 
models  of  communication  channels  in  urban  communities,  learning  about 
convolutional  codes,  their  implementation  and  performance,  and 
deciphering  the  Viterbi  algorithm  and  the  associated  bit  error  rate 
bounds.  This  research  continued  throughout  every  aspect  of  this 
thesis.  The  list  of  references  in  the  Bibliography  contain  very 
good  sources  for  topics  in  this  area  of  communications  and  other 
related  areas. 

This  report  was  limited  in  scope  to  two  specific  encoders  and 
four  specific  jammers  for  ease  of  model  development  and  implementation 
in  a  simulation  program.  However,  the  intent  was  not  to  restrict 
the  initial  problem  to  just  these  components.  The  ability  tc  add 
still  another  jammer  or  even  the  encoder-decoder  operation  exists. 

The  simulation  program  reflects  the  characteristics  and  statistics  of 
an  additive  white  Gaussian  noise  channel. 
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Electronic  countermeasures  (ECM) ,  such  as  broadband,  switched 
broadband,  continuous  wave  and  multitone  jamming,  are  designed  to 
disrupt  communication  systems.  The  subject  of  this  thesis  is  the 
effect  of  these  ECM  on  the  performance  of  error  correcting  codes  in 
an  additive  Gaussian  noise  channel.  The  channel  is  modeled  with  two 
convolutional  coders  having  a  rate  of  \  and  constraint  lengths  of  2 
and  6,  an  8-ary  FSK  modulator  and  demodulator,  and  a  decoder  based 
upon  the  Viterbi  algorithm.  Two  measures  of  effectiveness  are  the 
probability  of  bit  errors  P^Ce)  for  the  waveform  channel  and  the  bit 
error  rate  (BER)  for  the  overall  channel  performance.  It  is  shown 
that  both  (g)  and  BER  are  increased  by  the  previously  mentioned 
ECM. 

To  aid  in  evaluating  the  performance  of  these  coders,  a 
simulation  program  was  written.  In  the  program,  the  transmission 
signal,  the  channel  noise,  and  the  jamming  signal  are  independent 
of  each  other.  The  basic  conclusions  for  these  two  convolutional 
codes  are  that  the  dual -three  encoder  with  the  longer  codeword 
performs  better  in  the  no  jamming  or  in  a  Gaussian  jamming  environment, 
while  the  CC1  encoder  with  the  shorter  codeword  performs  better  in  a 
CW  jamming  environment. 


I .  INTRODUCTION 


Background 

Digital  communications  have  grown  more  popular  in  recent  years 
for  many  reasons.  One  such  reason  is  that  digital  transmissions ,  as  a 
subset  of  digital  communications,  provides  a  greater  flexibility  that 
is  not  available  with  analog  transmission.  Certain  advantages  of 
digital  over  analog  transmission  are  1)  for  analog  transmission  the 
receiver  attempts  to  trade  the  original  waveform  with  as  high  of 
fidelity  as  possible;  whereas,  for  digital  signaling  the  receiver 
decides  which  of  a  set  of  finite  signals  was  sent,  making  the 
probability  of  error  the  measure  of  effectiveness;  2)  digital 
regenerative  repeaters  reconstruct  the  message  signal;  whereas, 
analog  repeaters  amplify  the  signal  and  noise;  3)  digital  represen¬ 
tation  is  a  more  flexible  form  since  all  digital  signals  can  be  handled 
the  same  way  in  a  channel;  and  4)  a  digital  transmission  system  uses 
digital  integrated  circuits,  which  have  become  cost  effective  in 
electronic  circuit  design  (19:703). 

Two  basic  features  of  a  digital  communication  system  are  its 
signal  detection  ability  and  its  estimation  of  signal  parameters.  The 
probability  of  error  in  these  estimated  parameters  is  an  adequate  measure 
of  the  system's  performance.  In  signal  detection,  the  receiver  reduces 
the  received  waveform  to  a  set  of  numbers,  each  statistically  independent 
of  the  others.  A  correlation  operation  physically  generates  this  set  of 
numbers  invariant  to  the  decision  criteria.  The  second  feature,  estimation 


of  the  signal  parameters,  transforms  this  set  of  numbers,  identifying 
a  point  in  the  decision  space,  into  an  approximation  of  the  transmitted 
information.  Therefore,  once  this  set  of  numbers  is  formed  the  actual 
received  waveform  is  no  longer  important. 

The  major  objective  of  any  digital  communication  system  is  to 
minimize  the  total  probability  of  error.  Receivers  are  designed  based 
upon  expected  noise  and  other  interferences  to  reduce  the  error  probability 
as  much  as  possible.  Suppose  after  the  receiver  operation  is  optimized, 
an  extraneous  interfering  noise  is  received.  No  longer  is  the  receiver 
optimized  with  minimum  error  probability.  One  major  source  of  unwanted 
signals  in  a  military  environment  would  be  jamming  transmitters.  In 
communications,  electronic  countermeasures  (ECM)  are  designed  to 
disrupt  the  flow  of  communications  from  one  point  to  another.  This 
disruption  may  be  in  the  form  of  broadband  or  narrowband  noise  (25:17). 

The  effect  this  ECM  noise  has  on  the  communication  channel  and 
its  operation  are  addressed  in  this  report.  The  following  paragraph 


identifies  this  problem  with  a  specific  channel  model  and  four  jamming 


environments . 


Problem  Statement 

This  thesis  deals  with  the  analysis  of  the  effects  that  ECM  has  on 
error  correcting  codes  optimized  for  a  non-ECM  environment.  The  two 
major  portions  of  this  problem  are  first,  to  construct  an  additive  white 
Gaussian  channel  model  and  identify  four  types  of  jamming,  and,  second, 
to  predict  the  performance  of  two  convolutional  coding  schemes  in  each 
of  these  jamming  environments.  The  four  types  of  jamming  are  identified 
as  broadband,  switched  broadband,  continuous  wave,  and  multitone 


jamming.  Bit  error  probability  versus  signal -to-noise  ratio  is  the 
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major  criteria  for  comparison  of  the  two  coding  schemes.  A  simulation 
program  was  written  as  an  aid  in  evaluating  the  performance  of  this 
channel  model. 


General  Approach  and  Assumptions 

The  investigation  of  this  problem  began  with  research  of  published 
literature  describing  noisy  channel  models  (13:13,  20:11,  and  21),  the 
development  of  both  convolutional  codes  and  the  Viterbi  decoding 
algorithm,  and  ECM  techniques.  Two  convolutional  codes  were  then 
selected,  one  for  simplicity  and  the  other  for  its  longer  coding 
constraint  size.  The  channel  model  was  derived  mathematically,  assuming 
that  the  channel  noise  was  an  additive  white  Gaussian  process  and  that 
this  noise,  the  transmitted  signal,  and  the  jamming  interference  process 
were  all  statistically  independent.  The  8-ary  FSK  modulation  scheme 
is  a  well  known  and  well  used  scheme  for  transmissions  of  source  rates 
less  than  2400  bits  per  second  (30:1291).  The  demodulator  used  the 
ideas  of  detection  and  estimation  of  transmitted  signals  to  reconstruct 
the  coded  data  stream  for  the  decoder. 

The  analysis  consisted  of  looking  at  the  effects  of  these  jamming 
signals  on  the  probability  of  error  and  the  coding  bit  error  probabilities. 
The  basic  parameter  directly  affected  was  the  signal-to-noise  ratio 
which  in  turn  affected  the  other  quantities. 

A  simulation  program  was  constructed  from  the  mathematical  derivations 
of  the  channel  model.  The  structured  programming  approach  used  in  its 
development  allows  for  the  use  of  various  encoders,  decoders,  and 
modulation  schemes  without  adversely  affecting  the  program.  For  simpli¬ 
city  and  with  no  loss  in  generality,  an  all- zeroes  input  to  the  encoder 
was  used.  A  bit  error  out  of  the  decoder  was  then  a  "1"  (28:239). 


Prssentat  Lon 


The  organization  of  this  report  follows  the  development  of  the 
channel  and  jamming  models  and  their  analysis.  In  Cahpter  II  the 
different  components  of  the  channel  model  are  described.  The 
jamming  characteristics  are  defined  in  Chapter  III,  and  Chapter  IV 
deals  with  the  channel  performance  analysis  identifying  the  effects  of 
the  jamming  on  channel  model  parameters.  Chapter  V  contains  a 
description  of  the  simulation  program  and  the  eight  subroutines.  The 
conclusion  and  recommendations  are  discussed  in  Chapter  VI. 

In  addition,  three  Appendices  support  the  material  discussed  in 
the  report  Chapters.  Appendix  A  contains  the  mathematical  derivations 
in  support  of  Chapter  III.  Appendix  B  accomplishes  the  same  task 
for  the  bit  error  probability  discussed  in  Chapter  IV.  Finally, 
Appendix  C  contains  a  complete  listing  of  the  simulation  program. 


II.  CHANNEL  MODEL 


Introduction 

As  an  overview  f  digital  communications.  Figure  1  shows  the 
major  components  and  their  relationships  with  each  other.  Beginning 
with  the  source,  it  provides  discrete  source  symbols,  whether  digital 
data  or  digitized  analog  data.  The  representation  of  these  symbols 
will  be  discussed  later.  The  encoder  takes  these  source  symbols  and 
generates  channel  symbols  in  a  particular  manner  to  minimize  the 
effects  of  channel  noise.  At  this  point  the  modulator  inputs  the 
channel  symbols  to  the  channel  according  to  a  predetermined  modulation 
scheme.  The  channel  is  the  object  through  which  it  is  possible  to 
transmit  each  signal  during  a  specific  time  interval.  The  receiver, 
on  the  far  end  of  the  channel,  has  three  components:  demodulator, 
decoder  and  the  destination.  The  purpose  of  the  receiver  is  to  recreate 
the  original  source  symbols.  First  the  demodulator  performs  the 
inverse  operation  of  the  modulator  to  arrive  at  an  estimate  of  the 
original  channel  symbols.  The  decoder  then  maps  these  channel  symbol 
estimates  into  source  symbol  estimates  in  such  a  way  as  to  minimize 
the  effects  of  the  channel  noise  (19:704). 

To  initially  characterize  the  channel  model  for  this  thesis,  the 
source  provides  a  string  of  binary  digits  (...  b^,  b^  j,  b^  •••) 
where  each  bit  is  a  discrete  symbol  from  the  set  {0,1}.  The  encoder 
transforms  these  bits  into  codewords  of  length  two  or  six  depending 
upon  the  encoder.  The  elements  of  these  codewords  are  again  from  the 
set  {0,1}.  The  modulator  is  an  8-ary  FSK  system  which  selects  one  of 
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FIGURE  1:  DIGITAL  COMMUNICATION  SYSTEM  BLOCK  DIAGRAM 


eight  signals  for  every  three  bits  it  receives  from  the  encoder.  It 
then  inputs  this  signal  to  the  noisy  channel.  The  noise  component  of 
this  channel  is  an  additive  white  Gaussian  noise  (AWGN)  random  signal. 
Together,  this  noise  component  and  the  transmitted  signal  form  the 
received  signal.  The  receiver  is  a  maximum  a  posteriori  (MAP)  receiver. 
Its  first  component  is  the  demodulator  which  based  upon  the  received 
signal  determines  which  signal  of  the  modulator's  signal  set  (S^(t)} 
was  sent.  Then  it  reconstructs  the  channel  symbols,  three  bits  at  a 
time.  The  decoder  uses  a  shortest  route  algorithm  by  Viterbi  and  maps 
the  channel  symbols  into  an  estimated  bit  string  for  the  destination. 

The  effectiveness  of  this  channel  model  is  to  determine  the  number  of 
discrepancies  there  are  between  the  estimated  bit  string  and  the  original 
bit  string  from  the  source. 

The  ensuing  paragraphs  discuss  this  channel  model  in  more  detail. 
Both  encoding  schemes  are  discussed  as  well  as  the  modulation-demodu¬ 
lation  process  and  the  noisy  Gaussian  channel.  To  begin,  the  two 
encoders  will  be  discussed  showing  similarities  and  differences. 

Encoder 

The  encoder  observes  the  source  bits  at  its  input  and  generates  a 
sequence  of  bits  or  codewords  at  its  output.  This  generation  is 
accomplished  either  by  some  form  of  block  coding  or  some  type  of 
convolutional  coding  scheme.  The  encoders  in  this  model  are  two  types 
of  convolutional  coders.  Basically,  convolutional  coding  attempts  to 
avoid  exponential  growth  in  the  decoder  complexity  as  the  length  of  the 
codeword  increases  (28:227).  Convolutional  codes  can  be  analyzed  by 
several  methods.  However,  the  approach  taken  here  will  look  at  the 
shift-register  representation,  the  state  table,  and  the  state  diagram 
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of  each  encoder.  The  shift-register  diagram  illustrates  how  the 
individual  bits  are  combined  to  arrive  at  a  codeword.  The  state  table 
is  a  functional  table  that  describes  the  operation  of  the  encoder. 

It  contains  the  complete  information  about  the  present  state,  the  input 
bits,  the  output  codeword,  and  the  next  state.  All  possible  inputs  and 
present  state  combinations  are  enumerated.  The  state  diagram  is  a 
graphical  presentation  of  the  data  in  the  state  table.  The  circles 
identify  each  state,  and  the  directed  lines  between  states  show  the 
possible  state  transitions.  Input  and  output  data  accompany  the 
transition  lines. 

The  first  of  these  two  encoders  is  known  as  a  convolutional  code  1 
(CC1)  because  of  its  simplicity  (11:200).  Figure  2a  shows  the  shift- 
register  diagram,  Figure  2b  contains  the  state  table,  and  Figure  2c 
illustrates  the  state  diagram.  In  the  shift-register  diagram  the  two 
right  most  bits,  bi  ^  and  b^  2,  are  the  memory  bits,  and  they  identify 
the  four  possible  states  of  the  encoder.  Each  new  bit,  b^,  contributes 
to  the  generation  of  a  codeword,  then  it  shifts  into  position,  i-1. 

The  bit,  b^_j,  is  shifted  at  the  same  time  into  position,  i-2.  Once 
an  information  bit  has  occupied  all  three  positions  and  has  contributed 
to  the  generation  of  three  codewords,  it  is  shifted  out  of  the  register 
and  is  not  used  again.  The  state  table  for  CC1  shows  the  inputs  and  the 
codewords  for  the  possible  state  transitions.  Note  that  only  half  of 
the  table  is  filled.  This  results  from  the  fact  that  only  one  bit  is 
shifted  into  the  encoder,  but  two  bits  identify  the  encoder  state.  In 
the  state  diagram,  the  directed  lines  between  states  show  the  possible 
transitions  of  the  state  table.  If  the  transition  line  is  dashed,  then 
the  input  bit,  b^,  is  a  "1";  otherwise,  it  is  a  "0".  The  output  codeword 
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(c^c.,)^  is  in  parentheses  along  each  transition  line.  Therefore,  this 
encoder  maps  each  new  bit  with  the  previous  two  bits  to  arrive  at  a 
two-bit  codeword. 

The  second  convolutional  encoder  is  known  as  a  dual -three  encoder 
(2:42).  Figures  3a,  b,  and  c  contain  its  shift-register  diagram,  state 
table,  and  state  diagram,  respectively.  Again,  in  the  shift-register 
diagram  the  right  most  three  bits,  b^  ^ ,  b^_2>  and  b^  are  the  memory 
bits.  They  also  represent  the  present  state  of  the  encoder.  Since 
there  are  three  bits,  this  encoder  has  eight  possible  states.  The  three 
input  bits,  b^+2>  ^i+i*  aiu*  b^,  contribute  to  the  generation  of  the 
codeword,  (c^c^c^c^c^c^J  .  Then  they  are  shifted  into  positions,  i-1, 
i-2,  and  i-3,  and  become  the  next  state  of  the  encoder.  Again,  the 
usefulness  of  these  bits,  b^  ^ ,  b^  and  b.  is  over  after  they 
contributed  to  the  generation  of  two  codewords.  They  are  shifted  out 
of  the  register  and  are  not  used  again.  Referring  to  the  state  table, 
notice  that  all  possible  transitions  occur.  This  is  a  result  of 
replacing  all  the  memory  bits  with  the  input  bit  after  forming  each 
codeword.  The  output  codewords  are  in  octal  representation.  The  state 
diagram  for  this  encoder  is  very  detailed.  Since  each  state  has  eight 
possible  transitions,  only  those  transitions  from  state  "Oil"  are  shown. 
However,  every  state  has  the  same  configuration. 

This  encoder  is  more  complicated  than  the  first  one,  CC1.  The 
three  input  bits  and  six-bit  codewords  create  all  64  possible  transitions 
whereas,  in  the  CC1  encoder  one  input  allowed  only  half  of  the  state  tabl 
to  be  filled.  To  sum  up  this  operation,  the  encoder  maps  a  discrete  set 
of  inputs  into  a  codeword  of  particular  length.  In  the  first  encoder, 
this  discrete  set  was  three  bits  and  the  codeword  length  was  two;  in  the 


second,  the  discrete  set  was  six  bits  and  the  codeword  had  six  bits. 
In  both  cases  the  codeword  is  sent  to  the  channel  modulator. 


Channel 

The  channel  can  be  thought  of  as  the  modulator,  the  transmission 

media,  and  the  demodulator.  This  channel  is  called  a  "discrete-time 

memoryless  additive  Gaussian  channel  with  an  average  power  constraint" 

(11:90).  Its  inputs  come  from  the  set  of  these  binary  digits,  as  do 

its  outputs.  The  noise  components  are  independent,  identically 

distributed  normal  random  variables  with  mean  0  and  variance  N  /2.  To 

o 

begin  a  discussion  of  the  channel,  the  modulator  is  described  and 
characterized. 


Modulator.  The  modulator  tor  this  model  is  an  8-ary  frequency- shift 
keying  (FSK)  modulator.  For  every  three  bits  it  accepts,  the  modulator 
generates  one  of  eight  orthogonal  signals  for  transmission. 

Associated  with  this  modulator  is  a  set  of  eight  orthonormal  signals, 
An  orthonormal  set  of  signals  over  a  time  interval  T  seconds 
long  is  defined  in  the  following  manner  (22:33): 


2-1 


where 

$i(t)  =  ^==Sin  (i  ^  t)  i  =  1,  2,  ...  7,  8 

<$„  is  the  delta  function,  having  two  values  as  defined  above. 

The  decision  to  pick  one  of  these  signals  is  based  upon  the  three  input 
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bits  from  the  encoder.  Three  bits  generate  eight  octal  digits,  thereby 

creating  a  one-to-one  correspondence  between  the  three  bits  and  the 

orthonormai  signal  set,  {5^(t)}?_^.  Following  terminology  in  detection 

theory,  each  three-bit  input  can  be  thought  of  as  a  particular  hypothesis, 

II. .  For  eight  possible  combinations  of  ”1"  and  "0"  in  three  bits,  one 

hypothesis  identified  each  three-bit  input.  The  subscript  would  be  the 

octal  representation  of  the  input,  i.e.,  if  the  input  were  "100",  then 

H.  would  be  the  hvnothesis. 

4 

Based  on  the  set,  {i\(t)}?_^,  the  three-bit  inputs,  and  knowing 
that  all  transmitted  signals  have  equal  energy,  E,  the  modulator  forms 
the  set  of  transmission  signals,  {S^ ( t ) } ®_ 1 .  These  signals  are: 

Si(t)  =  /F  Sin  (i  p-  t)  i  =  l,  2,  ...  8  2-2 

where 


t  belongs  to  [0,T] 

The  modulator's  work  is  complete  when  it  forms  a  transmission 
signal  from  a  three-bit  input.  As  this  signal  is  transmitted  into  the 
noisy  channel,  the  modulator  accepts  three  more  bits  and  begins  forming 
another  transmission  signal.  Every  T  seconds  an  S^(t)  is  transmitted 
into  the  channel.  The  noise  in  this  channel  is  the  subject  of  the  next 
paragraph . 

Gaussian  Noise  Source.  This  noise  source  is  a  "discrete-time ,  memoryless 

Gaussian  source"  where  the  source  alphabet  is  the  set  of  real  numbers, 

and  the  output  is  a  sequence  of  independent  identically  distributed 

normal  random  variables  with  mean  0  and  vaiance  S,N'  (11:94).  The  N 

o  o 
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term  is  the  noise  spectral  density.  Two  reasons  for  selecting  an  AWGN 
source  are:  1)  it  has  minimum  complexity  and  2)  it  is  an  accurate  model 
for  an  important  class  of  communication  systems  (28:50). 

The  model  is  depicted  by  a  summing  junction  (Figure  1)  where  n(t) 
is  a  stationary  random  Gaussian  process  whose  power  density  spectrum  is 
constant  over  a  bandwidth  larger  than  the  front-end  bandwidth  of  the 
receiver.  To  form  the  received  signal,  r(t) ,  the  noise  process,  n(t) 
is  added  to  the  transmitted  signal,  S^(t). 

r(t)  =  n(t)  +  S. (t)  i  =  1,  2,  ...  8  2-3 

The  received  signal,  r(t) ,  is  the  input  to  the  receiver's  demodulator. 

Demodulator.  Basically,  the  demodulator  performs  the  inverse  operation 
of  the  modulator  (see  Figure  4) .  By  taking  r(t)  and  decomposing  it  over 
the  interval  [0, T]  with  the  set  of  orthonormal  functions,  <J>m(t),  the 
r^  coefficients  become  statistically  independent  Gaussian  random 
variables  with  equal  variances,  4N  (26:197).  From  this  coefficient, 
r^,  the  demodulator  attempts  to  estimate  the  correct  hypothesis,  1L, 
in  the  presence  of  interfering  noise.  Then  from  this  hypothesis,  an 
estimate  of  the  encoded  bit  string  is  formed,  three  bits  at  a  time. 
Therefore,  the  demodulator  has  two  main  operations.  The  first  is  to 
decompose  r(t)  into  eight  random  coefficients  by  using  the  orthonormal 
functions  in  {$u(t)}?_j,  and  the  second  operation  is  to  construct  the 
estimate  of  the  encoded  bit  string. 

To  accomplish  this  first  operation,  recall  Equation  2-3;  the 
received  signal,  r(t),  is  the  sum  of  the  transmitted  signal  and  the 
channel  AWGN  signal.  The  transmitted  signal  set  was  formed  from  the 
product  of  an  orthonormal  set  of  functions  (4u(t)}?_j  and  the  signal 
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energy.  This  set  {$^(t)}?_j  is  a  system  of  basis  functions  which  with 
the  property  of  Equation  2-1  span  the  entire  range  of  values  of  r(t). 
Therefore,  this  set  of  eight  equally  likely,  equal-energy  and  mutually 
orthogonal  functions,  {S^(t)}?_j,  partitions  the  range  space  of  r(t) 
into  eight  decision  regions.  Each  decision  region  in  the  partitioning 
has  a  one-to-one  correspondence  with  the  eight  message  signals;  such 
that,  if  r(t)  is  determined  to  be  in  one  of  these  regions,  then  the 
corresponding  S^(t)  is  assumed  to  be  the  transmitted  signal. 

To  find  where  r(t)  is  in  its  range  space,  eight  coefficients  are 
calculated.  These  coefficients  are  like  weighted  values  that  indicate 
how  well  r(t)  correlates  with  each  of  the  orthonormal  functions,  $^(t), 
over  the  time  interval  [0,T] .  The  largest  coefficient  identifies  a 
particular  basis  function,  <J>^(t),  which  in  turn  corresponds  to  a  decision 
region  and  a  particular  transmitted  signal,  S^(t).  These  eight  co¬ 
efficients  are  formed  in  the  following  manner: 


where 


r^  is  the  received  random  variable  in  the  mth  region 

n  is  the  AWGN  variable 

m 

6.  ,  is  defined  in  Equation  2-1 

l ,  m+  i 
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( 


E  is  the  signal  energy 

The  rm  coefficient  is  now  the  sum  of  two  components,  each  of  which 

becomes  a  weighted  value  in  itself  of  the  correlation  between  4>m+^(t) 

and  the  component  signals,  i.e.,  S^(t)  and  n(t).  Since  S^t)  and  n(t) 

are  independent  of  each  other  and  n(t)  is  an  AWGN  random  signal,  the 

r^  coefficients  are  statistically  independent  Gaussian  random  variables 

with  a  variance  equal  to  I5N  .  The  mean  of  r  depends  upon  6. 

n  o  m  r  r  i,m+l 

Equation  2-S  summarizes  the  probability  density  function  of  r^  given  a 
particular  signal  was  transmitted,  S^t),  or  a  particular  hypothesis, 

H.,  is  true. 

N(v^,  %N  )  m=i 

0  2-5 

N (  0,  m/i 

To  finish  this  first  operation,  the  demodulator  selects  the  largest  r 

m 

coefficient. 

The  second  operation  is  to  recover  the  estimated  coded  bit  string 

from  the  largest  r  coefficient.  This  specific  r  identifies  the 
m  r  m 

particular  $m+1(t)  that  correlated  the  most  with  r(t)  .  From  $>m+^(t), 
the  decision  region,  the  assumed  transmitted  signal,  Sm+j(t),  and  the 
particular  hypothesis,  H  ,  are  determined.  Since  each  hypothesis 
identifies  one  octal  digit,  the  demodulator  can  now  estimate  the  three- 
bit  input  to  the  modulator.  These  bits  are  held  in  a  buffer  until  a 
complete  code'  ’  can  be  passed  on  to  the  decoder.  This  buffer  is 
necessary  since  codewords  are  fractionally  prepared  for  transmission  in 
the  modulator. 

Probability  of  Error.  The  decision  of  choosing  an  hypothesis  based  upon 


fr  ,H,W 
m  1 
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r  will  lead  to  a  correct  decision  most  of  the  time.  Two  measures  of 
m 

decision  error  and  correctness  are  their  respective  probabilities,  P(e) 

and  P(C).  In  choosing  the  largest  r^,  the  receiver  compares  two 

coefficients  and  retains  the  larger  one.  This  pairwise  comparison  is 

done  until  the  largest  r  is  found.  Rather  than  trying  to  visualize  an 

m 

eight  dimensional  space,  one  comparison  will  be  analyzed.  Figure  5 

shows  the  r^  axis  with  the  two  probability  density  functions  (Equation  2-5) 

for  r^  and  r^,  given  is  true.  Let  r^  equal  a  (see  Figure  5).  Then, 

for  a  correct  decision,  r^  must  be  in  the  interval  I  ,  meaning  r^<a, 

and  for  an  incorrect  decision,  r.  must  be  in  the  interval  I  which  means 

1  w 

r^>a.  Equation  2-6a  follows: 

P(C|Hk)  =  P(rielc|Hk)  =  P(ri<a|Hk)  2-6a 

=  Jl  fr.|H„Wdx 

c  i 1  k 

Likewise,  the  P(e)  follows  from  the  fact  it  equals  1  -  P(C). 

P(e|Hk)  =  1  -  Jj  fr  iH  (x)  dx  =J:  ff  ,H  (x)  dx  2-6b 


c  i 1  k 


w  i 1  k 


=  P(r.>a|Hk) 


Moving  from  the  binary  comparison  case  to  the  problem  of  determining  the 

P(C)  and  P(e),  the  P(CjHk,rk)  for  the  8-ary  basis  system  means  that  every 

coefficient  is  less  than  rk  or  P(rQSrk,  r^r^  ...  r7srk) .  Since  the 

r  's  are  statistically  independent  and  identically  distributed,  this  last 
m 

probability  becomes  the  product  of  the  probabilities  of  the  individual 
inequalities,  [P(r^<rk)]7.  For  eight  equally  likely,  equal-energy  signals, 
by  multiplying  by  the  probability  distribution  function  of  r^  and 
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FIGURE  5:  DECISION  SPACE  FOR  THE  DEMODULATOR 


integrating  over  the  range  of  r  ,  the  P(C|H^'!  is  calculated.  Van  Trees 
(26:261-3)  has  the  details  for  finding  these  probabilities.  The  results 
for  P(C)  and  P(e)  are  Equations  2-7a  and  2-7b: 


P(C)  =•  (2tt)  _ls  j  exp[-^(x->  2E/Nq)  2]  [  (2tt)  j  e  ‘^da^dx  2-7a 

—  CO  -00 

P(e)  =  1  -  P(C)  2-7b 

where 

Tk 

x  =  - 

Ai  /  2 
o 

Returning  to  the  explanation  of  the  components  of  Figure  1,  the 
decoder  receives  the  estimate  of  the  codeword  from  the  demodulator. 

What  happens  to  this  codeword  is  discussed  in  the  description  of  the 
decoder. 

Decoder 

The  basic  function  of  the  decoder  is  to  reverse  the  mapping  of  the 
coded  message  into  the  original  bit  string  with  minimum  error.  The  use 
of  the  Viterbi  decoding  algorithm  in  convolutional  coding  appears  to  be 
the  most  cost-effective  forward  error  control  technique  for  medium 
accuracy  requirements.  The  algorithm  is  relatively  simple  to  implement 
and  is  effective  since  it  essentially  implements  maximum- likelihood 
decoding  (19:716).  Maximum- likelihood  decoding  means  that  the  decoder 
compares  the  received  codeword  with  all  possible  encoder  state 
transition  outputs  and  selects  the  transition  whose  coded  output  most 
correlates  with  the  received  codeword.  This  correlation  is  based  upon 
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a  bit-by-bit  comparison  counting  the  number  of  bits  that  differ  in 
both  the  received  codeword  and  each  transition  output.  The  number  of 
bits  that  differ  is  referred  to  as  the  Hamming  distance  between  two 
binary  sequences. 

Basically,  this  decoder  attempts  to  simulate,  with  the  received 
codeword,  the  behavior  of  the  encoder.  With  the  help  of  the  enccJ  r's 
state  diagram,  the  encoder's  behavior  can  be  traced  from  state  to  state 
traveling  along  the  transition  lines  associated  with  the  minimum 
Hamming  distances  relative  to  the  received  codeword.  Difficulty  in 
keeping  track  of  each  state  transition  occurs  when  the  same  transition 
is  used  several  times.  A  trellis  diagram  is  a  modified  state  diagram 
that  has  included  a  dimension  of  time  (11:207).  Figure  6  illustrates 
the  construction  of  a  general  trellis  diagram.  Each  row  represents 
each  state  in  the  state  diagram,  and  the  transition  paths  correspond 
to  the  same  transition  lines  in  the  state  diagram.  The  input  bits  and 
the  output  codeword  accompany  each  transition  path.  Another  version  of 
the  trellis  which  is  used  later,  contains,  instead  of  each  output  code¬ 
word  along  the  transition  paths,  the  Hamming  distance  between  the 
received  codeword  and  the  transition  output  codeword.  The  transition 
depth  is  the  added  time  dimension  making  this  diagram  different  from 
the  state  diagram. 

To  traverse  the  shortest  cumulative  Hamming  distance  through  the 
trellis,  the  following  algorithm  is  used  (11:211): 

First,  the  three  basic  variables  are  defined: 

Uj (a)  =  the  cumulative  Hamming  distance  at  depth  j  and  state  a 
B  (a)  =  the  cumulative  output  bit  string  at  depth  j  and 
state  a 
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of  the  information  bit  string.  Therefore,  upon  reconstructing  the 
information  bit  string,  the  last  two  bits  can  be  discarded  since  they 
were  not  part  of  the  original  information  bit  string. 

Referring  to  the  CC1  trellis  diagram  in  Figure  7a,  a  straight 
forward  way  to  decode  codewords  would  be  to  calculate  the  Hamming  distances 
for  the  received  codeword  on  all  paths  through  the  trellis.  The  path 
with  the  smallest  distance  would  then  be  selected,  and  the  information 


bits  corresponding  to  the  path  would  form  the  decoder  output.  However, 
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the  number  of  paths  for  an  L  bit  information  sequence  is  2^.  For 
large  L  this  approach  becomes  impractical  (9:839). 

The  Viterbi  algorithm  greatly  reduces  the  number  of  paths  to  check 
by  taking  advantage  of  the  special  structure  of  the  trellis.  Figure  7a 
shows  a  definite  periodic  form  after  reaching  depth  3.  After  this 
depth,  each  of  the  four  states  are  entered  from  only  two  preceding 
states.  The  decoder  then  calculates  the  Hamming  distance  associated  with 
each  of  the  two  paths  entering  a  given  state  and  eliminates  the  path 
that  has  the  greater  distance  from  further  consideration.  This  decision 
process  is  done  four  times  during  each  trellis  depth.  After  the  decoding 
process  is  complete,  only  one  path  remains  leading  to  each  state.  The 
decoder  then  proceeds  to  the  next  trellis  depth  and  repeats  this  process. 

For  the  trellis  in  Figure  7a,  there  are  eight  paths  at  depth  3. 
Decoding  at  depth  3  eliminates  one  path  entering  each  state.  The  result 
is  that  four  paths  are  left.  Proceeding  on  to  depth  4,  again  the  decoder 
has  to  reduce  eight  paths  to  four. 

Now,  the  assumption  of  finishing  in  the  "00"  state  aids  in  selecting 
a  single  most  likely  path  having  the  shortest  Hamming  distance.  Without 
this  assumption  the  decoder  could  never  settle  on  one  most  likely  path. 

It  always  retains  one  path  to  each  state  after  the  decoding  process. 

There  would  be  four  more  likely  paths  at  any  depth  greater  than  3. 

Forcing  the  encoder  to  the  "00"  state  by  inputting  two  zeroes  to  the 
encoder  after  th  source  input  string  is  the  way  this  decoder  can  select 
a  single  most  likely  path  having  the  smallest  Hamming  distance.  This 
optimum  path  finishes  in  the  prearranged  state  of  "00"  (9:839). 

Figure  7b  shows  the  four  more  likely  intermediate  paths  and  the  single 
minimum  distance  path  from  depth  0  to  depth  6.  The  information  bit  string 
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for  the  optimum  path  is  written  below  the  trellis. 

Applying  this  basic  algorithm  to  the  second  decoder,  this  decoder 
assumes  again  that  the  encoder  began  and  ended  in  the  "00"  state. 
Referring  to  the  dual-three  trellis  diagram  in  Figure  8a,  a  straight 
forward  way  to  decode  the  received  codewords  would  be  very  involved 
since  between  each  trellis  depth  there  are  64  transition  paths  after 
depth  2.  The  number  of  paths  grows  exponentially  from  one  depth  to 
another.  Here  the  Viterbi  algorithm  can  greatly  reduce  the  workload  of 
the  straight  forward  method.  After  depth  2  the  trellis  shows  a 
definite  periodic  pattern.  Each  of  the  eight  states  are  entered  from 
the  eight  preceding  states.  The  decoder  then  calculates  the  Hamming 
distance  associated  with  each  of  the  eight  paths  entering  any  given 
state  and  retains  only  the  one  path  with  the  least  Hamming  distance. 

The  others  are  discarded.  This  seven-path  elimination  is  done  eight 
times,  once  for  each  state,  during  each  depth  of  the  trellis.  After 
the  decoding  process  is  complete,  only  one  path  remains  leading  to 
each  state.  The  decoder  then  proceeds  to  the  next  trellis  depth  and 
repeats  the  process,  reducing  64  paths  to  8. 

Now  the  assumption  of  finishing  in  the  "000"  state  aids  in 
selecting  a  most  likely  path  having  the  shortest  Hamming  distance. 
Without  the  prearranged  finishing  state,  the  decoder  would  always  retain 
eight  paths  from  one  depth  to  another.  Figure  8b  shows  the  eight  more 
likely  intermediate  paths  and  the  single  minimum  distance  path  from 
depth  0  to  depth  4  (9:839).  The  information  bit  string  for  the  shortest 
path  is  written  below  the  trellis. 

Destinat ion 

The  final  portion  of  Figure  1,  the  destination,  receives  the 
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Received 

Codewords 


001101 


011011 


101001 


000100 


State 


8b.  OPTIMUM  DECODING  PATH 

FIGURE  8:  TRELLIS  DIAGRAMS  FOR  DUAL-THREE 
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reconstructed  bit  string.  This  bit  string  is  the  estimation  of  the 
original  input  bit  string  to  the  encoder.  The  receiver  is  constructed 
with  a  minimum  probability  error  criterion  and  equal  a  priori 
probabilities  of  the  hypotheses  (26:257).  Since  all  eight  signals 
have  equal  energy,  the  receiver  simply  chooses  and  processes  the  signal 
that  is  the  most  correlated  with  r(t) .  Processing  and  decoding 
translate  the  most  correlated  signal  into  the  estimate  of  the  original 
bit  string. 


*  • 
f 
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III.  JAMMING  MOGULS  AND  PERFORMANCE  ANALYSIS 


The  strategy  of  a  jammer  is  to  confuse  the  receiver  and  make  it 
incorrectly  determine  the  transmitted  message.  In  creating  this  error, 
the  jammer  operates  at  the  frequencies  of  the  communication  system. 

This  Chapter  discusses  four  types  of  jamming:  broadband,  switched 
broadband,  continuous  wave,  and  multitone  jamming.  The  following  models 
consider  all  jamming  signals  as  an  additive  interference.  Figure  9 
shows  the  additive  jamming  signal  in  the  Digital  Communication  System 
of  Figure  1.  The  new  parameter  in  this  Figure  is  J(t),  the  interfering 
noise  process  which  is  independent  of  both  the  signal  and  n(t).  Depending 
upon  the  type  of  jamming,  J(t)  is  either  deterministic  or  random.  This 
interference  process  is  further  characterized  in  the  following  sections. 

Broadband  Jamming 

Broadband  jamming  is  the  simplest  jamming  technique  to  envision 
when  it  has  a  uniform  spectral  density  over  a  bandwidth  greater  than  or 
equal  to  the  bandwidth  of  the  receiver.  The  effect  in  the  receiver  is 
another  additive  white  Gaussian  noise,  which  is  a  zero-mean  random 
process.  In  Equation  3-1,  the  received  signal  is  the  sum  of  the  transmitted 
signal,  the  channel  AWGN  signal  and  the  jamming  signal. 

r(t)  =  /ES^t)  +  n  (t)  +  J(t)  te[0,T]  3-1 

where 

E  is  the  received  signal  energy 
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FIGURE  9:  DIGITAL  COMMUNICATION  SYSTEM  WITH  ADDITIVE  INTERFERENCE,  J(t) 


S^(t)  is  the  message  signal  for  hypothesis  ^ 
n(t)  is  the  channel  AWGN  process 
J(t)  is  the  random  Gaussian  process 
T  is  the  signaling  time  interval 

From  this  set  of  message  signals,  {S  (t)},  and  from  the  fact  that  each 
signal  is  equally  likely  and  has  equal  energy,  the  minimum  attainable 
error  probability,  P(e),  is  easy  to  calculate.  This  set  of  message 
signals  is  a  particular  set  of  eight  orthogonal  signals,  or  a  system 
of  basis  functions,  having  the  property  of  Equation  2-1  over  the 
interval  [0,T]  and  spanning  the  entire  range  of  r(t)  .  With  these 
eight  equally  likely,  equal -energy  and  mutually  orthogonal  functions, 
a  decision  region  for  each  can  be  formed.  Each  decision  region 
identifies  a  partition  of  the  range  of  values  for  the  received  signal, 
r(t) .  Also,  these  regions  have  a  one-to-one  correspondence  with  the 
eight  message  signals;  such  that,  if  r(t)  is  determined  to  be  in  one 
of  these  regions,  then  the  corresponding  S^(t)  is  assumed  to  be  the 
signal  transmitted. 

To  determine  where  r(t)  happens  to  be  in  its  range  space,  eight 
coefficients  are  calculated.  Each  coefficient  is  like  a  weighted  value 
that  indicates  how  well  it  correlates  with  each  of  the  eight  basis 
functions  over  the  time  interval  [0,T].  The  largest  coefficient 
identifies  a  particular  basis  function  which  in  turn  corresponds  to 
a  decision  region  and  thus  a  particular  signal,  S.(t).  These  eight 

i 

coefficients  are  formed  in  the  following  manner: 
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r  is  the  received  random  variable  in  the  mth  region 
m  * 

n  is  the  AWGN  random  variable 
m 

j  is  the  jamming  random  variable 

p.  ,  is  the  correlation  coefficient  between  S.(t)  and  S  .  (t) 
i.nw-1  i  m+1 

The  r^  coefficient  is  a  sum  of  three  components,  each  of  which  becomes  a 
weighted  value  of  the  correlation  between  the  component  signal  (i.e., 

S.  (t) ,  n(t)  or  J(t))and  the  mth  basis  function.  Since  S  ^ ( t ) ,  n(t), 
and  J(t)  are  independent  of  each  other  and  both  n(t)  and  J(t)  are  AIVGN 
random  signals,  the  rm  coefficients  are  statistically  independent 
Gaussian  random  variables  with  a  variance  equal  to  the  sum  of  the 
individual  variances:  (Nq  +  Jq)/2.  The  mean  of  each  r^  depends  upon 
which  message  signal  or  hypothesis  is  true.  That  is,  if  r^  is  calculated 
and  S^(t)  is  sent,  then  when  m=i  the  mean  is  SE  and  when  mj*i  the  mean 
is  0.  Equation  3-3  summarizes  the  probability  density  function  of  r 

m 

given  a  particular  hypothesis,  Ih  ,  is  true. 

N(v1f,  (N  +  J  )/2)  m=i 

f  |h  Cx)  =  0  °  3-3 

m1  i  N (  0,  (N  +  J  )/2)  m^i 

oo 

With  the  receiver  design  described  in  Chapter  II,  the  decision  of 
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choosing  an  hypothesis  based  upon  the  ’•eceived  signal  coefficients,  r^ 

for  all  m=(l,  1,  ...  7,  will  lead  to  a  correct  decision  some  of  the  time. 

This  amount  of  time  for  being  correct  can  be  descrioed  by  the  probability 

ot  correctness,  P(C),  and  vice  versa  for  the  probability  of  error,  P(e). 

In  choosing  the  largest  r  ,  the  receiver  compares  two  coefficients  and 
'  m 

retains  the  larger  one,  only  to  be  compared  with  the  next  coefficient 
until  this  pairwise  comparison  process  is  complete.  Rather  than  trying 
to  visualize  an  eight  dimensional  space,  one  comparison  will  be  analyzed. 
Figure  10  shows  the  received  coefficient  axis  with  the  two  probability 
density  functions  for  r ^  and  r^  given  hypothesis,  (1^,  is  true.  Let 
r^  equal  a  (see  Figure  10).  Then,  for  a  correct  decision,  r^  must  be 
in  the  interval  1  ,  meaning  r.<a,  and  for  an  incorrect  decision,  r. 
must  be  in  the  interval  which  represents  those  values  of  r^  greater 
than  r^.  Equation  3-4a  follows: 

P(C|Hk)  =  P(r.eIc|Hk)  =  F(r.<a|Hk)  3-4a 

=  /l  fr.K(x)dx 
c  i 1  k 

Likewise,  the  P(e)  follows  from  the  fact  it  equals  l-P(C). 

P (e)  =  1  -  Jj  fr  |H  (x)  dx  =  Jj  fr  ,H  (x)  dx  3-4b 

c  i '  k  w  i 1  k 


*  P(r^>a|Hk) 

Moving  from  the  binary  comparison  case  to  the  problem  of  determining 
the  P(C)  and  P(e),  the  P (C 1 1  ^ , rk )  for  the  8-ary  basis  system  means 
that  every  coefficient  is  less  than  rk  or  P(rQ<rk>  rj$rj(>  •••>  r^<r  ) . 
Since  the  r  's  ate  statistically  independent  and  identically  distributed, 
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FIGURE  10:  DECISION  SPACE  FOR  BROADBAND  JAMMING 


this  last  probability  becomes  the  product  of  the  individual  inequalities. 


[P(r.<r ,) ] 7 •  For  eight  equally  likely,  equal-energy  signals,  by 

1  K 

multiplying  by  the  probability  distribution  function  of  r  and  integrating 
over  the  range  of  r^  the  P(C|ll^)  is  calculated.  Appendix  A  has  the 
details  for  finding  this  probability.  The  results  for  P(C)  and  P(e) 
are  in  Equations  3-5a  and  3-5b: 


1  °°  [(X-/2E/ 

P(C)  =  /exp - 


i' wiu- 

2  JL»/2Tr-<» 


P(e)  =  1  -  P(C) 


3-5b 


where 


x  = 

Switched  Broadband  Jamming 

This  type  of  jamming  is  defined  as  the  product  of  a  periodic 
switching  function,  z(t),  and  the  broadband  jamming  signal,  J(t)  (see 
Figure  11).  Again,  Equation  3-2  applies;  however,  before  evaluating 
the  third  integral, 

T 

z(t)  J(t)  SCt)  dt, 

0  m 

statistical  information  about  z(t)  is  required.  The  average  duty  cycle, 
d,  is  defined  to  maintain  the  average  jamming  power  at  Jq/2  for  any 
time  interval.  The  probability  distribution  function  for  z(t)  is 
P  (z)  =  d  when  z(t)  =  1  and  P„(z)  =  1-d  when  z(t)  =  0.  Then,  E{z(t)}  =  d 
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and  the  variance  of  z(tj  is  d(l-dj.  Letting  equal  the  above 
integral,  its  expected  value  and  correlation  function  become 


Hi  }  =  0 

•  m 


3-6 


K . (m,k)  = 


'■jJ  d(d  - 


1 


Sin(  fm+l]4”nd)  )  m=k 


4^ (m+  1 ) n 


liJ  d( 


1 


2Tn(m-k) 


1 


27rn(m+k  +  2) 


Sin(  (m-k)2'rnd) 


Sin(  (m+k  +  2)  27tnd)  )  m^k 


Appendix  A  contains  the  mathematical  derivation  of  these  parameters. 

Figure  11  illustrates  the  subsystem  that  produces  the  interference. 
z(t)J(t).  Both  functions,  z(t)  and  J(t),  are  independent  random 
processes.  This  operation  is  a  linear  operation  on  a  Gaussian  random 
process,  J(t);  therefore,  the  result  is  a  Gaussian  random  process. 


z(t) 


J(t)  A  J(t)z(0 

- (x) - 


FIGURE  11:  JAMMING  SUBSYSTEM 


Letting  j  equal  the  third  integral  in  Equation  3-2  again,  the 

distribution  of  the  coefficients,  r  ,  is  similar  to  the  broadband  case 

m 

where  the  covariance  function,  K.(m,k),  of  Equation  5-7  is  used. 


fr  = 

m  i 


N (  0, 

l$N  ♦ 

K.lm.k)  4 

m/k 

0 

J 

N(  A, 

liN  + 
o 

k  j  (.  m ,  k )  ) 

m=k 

3-8 
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To  determine  the  probability  of  error,  P(e),  the  same  procedure 
is  used  as  for  the  broadband  case.  This  is  valid  since  the  total 
interference  is  a  Gaussian  process.  The  only  difference  between  the 
broadband  and  this  switched  broadband  case  is  the  covariance  function, 
K  (m,k).  In  the  latter  case,  the  P(e)  and  P(C)  are  as  follows: 


P(C) 


oo 

1  f 

- /  exp 

(x-»/E/ g2)2 

i  x 

- -  exp  (-'-ja2)  da 

-OO 

2 

/2tF 

-CO 

dx 


3-9 


P (e)  =  1  -  P(C) 


where 


r 

m 


g 2  =  '-sN  +  K  (m,  k)  ,  x  =  — rr 

°  J  /g2 


Continuous  Wave  Jamming 

This  jamming  is  a  signal  at  a  specific  frequency  and  phase  angle, 
i  .  e . : 

J(t)  =  OsJoK  Sin  (w.t  ♦  0)  3-11 

where 


J  is  the  jamming  power  density  parameter 
o 

h).  is  the  jamming  frequency 
6  is  the  phase  angle 

Referring  back  to  liquation  3-2,  the  third  integral  is  equal  to  j  . 
Depending  upon  whether  or  not  u>.  equals  the  frequency  of  S^+j(t), 
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h  vj  T  cos  0 
o 


idj  =  m2im/T 


/J~ 

2co 

0 

m 

T 

2  2 

a).  -  U) 

L  J  m 

cos  9  sin  (w. t) 


3-12 


+  sin  0  (cos  ukT  -  1) 


ojj  f  m27rn/T 


The  mathematical  derivation  from  j  is  in  Appendix  A.  The  phase  angle, 
9,  is  a  random  variable  assumed  to  be  uniformly  distributed  over  the 
interval  [--i,tt].  With  its  mean  equal  to  0  and  its  variance  equal  to 
n2/3,  the  expected  value  of  j  is: 


»s(j  T)  CO.  =  UL 

r .  i  o  i  m 

elj  J  =  /  >-  \ 

(hU/T)2  (2u  /  (cj- 2  -  co2))  sin  w.t)u.  * 
l\o  m  j  m  j  /  j  m 


3-13 


The  random  variable,  j  ,  becomes  deterministic  once  the  phase  angle  is 
known.  Thus,  jm  becomes  a  parameter  which  affects  the  mean  value  of 
the  received  signal,  r^.  The  probability  distribution  function  for  r 
is  N(/£T  m  +  Jm  >  Nq/2)  .  Figure  12  shows  the  effect  this  jamming  has 
on  the  decision  space  to  determine  a  correct  or  an  incorrect  received 
signal . 

The  calculations  for  the  probability  of  error  and  correctness  are 
straight  forward  and  can  be  followed  in  Appendix  A. 


P(C)  =  — f  exp 
/2tT  J 

-CO 

P(e)  =  1  -  P(C) 


/E  +  j 

’  X 

-*a(x  -  — - -)2 

f  — —  exp(-^a2)da 

/  2 

„  0  J 

J  /2m 

L  -oo  J 

dx  3-14 


3-15 


Multitone  Jamming 

This  jamming  technique  employs  more  than  one  continuous  wave 
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FOR  CW  JAMMING 


is  the  average  of  the 


jamming  signal.  The  jamming  random  variable,  j 
individual  signal  contributions  to  the  interference.  Therefore,  this 
variable  becomes: 

N  T  N 

j  =  ~  Z  /  J .  ft)  S  ,  (t)  dt  =  i-  z  v  .  3-16 

m  N  i  =  i  J0  1  m+l  N  i  =  1  m,i 

*c 

where 

N  is  the  number  of  tones 

J^(t)  is  the  ith  jamming  tone 

S  , (t)  is  the  m+lst  orthogonal  function 
m+ 1  & 

v  .  is  the  ith  jamming  tone  random  variable 
m,  l  .6 

Each  vm  ^  is  defined  as  the  continuous  wave  jamming  random  variable, 

jm»  in  Equation  3- Id.  As  in  the  continuous  wave  case,  j^  is  a  parameter 

that  affects  or.’.,  the  mean  of  the  received  coefficient,  r  ,  once  the 

m 

phase  is  known  Therefore,  the  probability  distribution  function  for 

r  is  N(vC  j  ♦  j  ,  N  / 2) ,  where  j  is  defined  in  Equation  3-16.  This 
m  .  - ■  Jm  o  m  M 

j  has  the  s;r:  characteristics  as  the  set  of  coefficients  of  Equation  3-1 
For  the  multi.*  no  case,  Figure  12  also  applies.  It  shows  the  effects 
multitone  jamming  has  on  the  decision  space  and  ultimately  on  the 
correctness  or  error  of  the  receiver  output. 

The  calculations  for  the  P(C)  and  P(e^  are  in  Appendix  A.  The 
results  are  the  same  as  in  Equations  3-14  and  3-15,  with  the  only 
exception  being  the  determination  o;  j  . 

Probability  of  Bit  Error 

In  the  communication  system  diagrammed  in  Figure  1,  the  P(e)  and 
P(C)  are  determined  from  the  received  signal  converted  to  a  3-bit  word 
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(i.e.,  000,  001,  ...,  111).  Thses  bits  are  assembled  for  the  decoder. 


However,  it  is  possible  to  have  one  or  two  of  the  three  bits  correct 
when  an  incorrect  decision  is  made  in  the  demodulator.  For  instance, 
assume  the  correct  receive  signal  was  "101".  Any  one  of  the  other 
seven  channel  symbols  is  a  channel  symbol  error.  Looking  at  only  the 
first  bit,  "1",  ii  the  incorrect  signal  were  "100",  "110",  or  "111", 
then  the  first  bit  would  be  correct.  Therefore,  the  P(e)  for  the 
received  signal  differs  from  the  probability  of  error  for  each  bit, 

P .(e) .  The  calculation  for  P  (e)  in  relation  to  P(e)  is  in  Appendix  A 
and  is  the  same  for  each  of  the  previously  discussed  jamming  techniques. 


Pb(e)  =  y  P (e)  3-17 

P. (C)  =1-4  P(e)  3-18 
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V.  CHANNEL  PERFORMANCE  .ANALYSIS 


Int roduct ion 

The  two  previous  Chapters  have  characterized  the  general  model 
of  the  digital  communication  system  in  Figure  1  and  have  identified 
four  types  of  jamming.  The  probabilities  of  symbol  and  bit  errors 
were  calculated  for  the  effects  of  each  jamming  type  on  the  recovery 
of  source  information  at  the  destination.  Specifically,  how  this 
jamming  affects  the  channel  model  is  the  subject  of  this  Chapter. 

First,  the  probability  of  bit  error  will  be  analyzed  in  relation  to 
the  signal -to-noise  ratio  as  it  is  affected  by  the  jamming,  and  then 
the  bit  error  probability  for  the  two  pairs  of  encoder-decoders  will 
be  bounded  tightly.  This  bound  will  reflect  the  effects  of  the 
jamming  on  the  coding  operation. 

Before  moving  to  the  probability  of  bit  error  analysis  some 
specific  characteristics  of  this  channel  are  required.  An  8-ary  FSK 
modulator-demodulator  design  is  the  basic  choice  of  designers  for  data 
speeds  up  to  2400  bits  per  second  (30:1291).  The  rate  for  data  transfer 
will  be  1800bps.  To  establish  a  bound  on  the  signal-to-noise  ratio  of 
the  channel,  signal  and  noise  thresholds  should  be  identified.  First 
of  all,  a  system's  noise  threshold  is  the  received  signal  level  at 
which  the  RMS  signal  is  equal  to  the  RMS  noise  at  the  input  to  the 
demodulator.  The  signal  threshold,  for  an  FM  system,  is  the  received 
signal  level  at  which  the  signal  peaks  begin  to  exceed  the  noise  peaks 
consistently.  At  this  point,  modulation  becomes  a  useful  tool  of 
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communications  (31:7).  The  FM  threshold  is  measured  10  dB  above  the 

noise  threshold.  For  this  channel  model’s  modulation  scheme,  the 

E 

signal -to-noise  ratio  (y^-)  will  be  greater  than  or  equal  to  10. 

‘  o 

Probability  of  Bit  Error  Analysis 

Figure  13  has  a  complete  listing  of  the  probabilities  of  bit 

error  for  each  jamming  case.  The  forms  are  slightly  different  from 

the  ones  in  the  previous  Chapters  to  identify  the  susceptibility  of 

P^Ce)  to  signal-to-noise  ratio  changes  (28:96).  Since  the  main  concern 

is  with  bit  errors,  several  of  this  channel's  parameters  will  refer 

to  "per  information  bit"  values.  Normalizing  the  probability  of 

symbol  error,  P(e),  to  its  "per  information  bit"  value,  P^fe)  equals 
4 

y  P(e).  Also,  since  Es  is  the  symbol  energy,  E^  becomes  its  "per 

information  bit"  value  with  Efa  equaling  y  Es  (19:711).  This  follows 

from  three  bits  generating  one  channel  signal  of  energy  E^  over  an 

interval  of  T  seconds.  For  the  no  jamming  situation.  Figure  14 

graphs  the  error  probability  for  eight  orthogonal  signals  (29:259). 

E  E 

Notice  that  the  abscissa  has  two  scales,  10  Log  s  and  10  Log?n  b. 

^^N-  luiT" 

The  relation  between  these  ratios  is  10  Log^  y  (-4.77  dB) .  Each 

still  identifies  the  same  points  on  the  two  curves,  P^(e)  and  P(e). 

For  the  no  jamming  case  the  noise  is  just  N  . 

Turning  to  the  broadband  jamming  case,  compare  the  two  equations 

in  Figure  13  for  no  jamming  and  for  broadband  jamming.  Both  equations 

have  the  same  basic  form,  but  the  signal  to  noise  is  smaller  in  the 

broadband  jamming  equation  since  the  noise,  N,  is  now  composed  of  both 

E 

the  channel  AWGN,  N  ,  and  the  jamming  AWGN,  J  .  Let  ^  equal  20 

J  E  o 

s 

(or  13  .  Now  as  J  increases,  ^  y  decreases,  and  the  probability 

o'o 
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of  error  increases.  Referring  to  Figure  14,  as  the  signal-to-noise 


ratio  decreases  the  P(e)  and  P^(e)  increase.  In  fact,  since  both 

noises  are  AWGN  processes,  the  effect  of  broadband  jamming  on  the  P^(e) 

of  this  channel  model  will  follow  the  curve  in  Figure  14.  However, 

E 

s 

when  J  equals  N  ,  10  Log,,,  —  is  10  dB,  which  is  the  modulation 
o  n  o  10  N 

threshold  discussed  earlier.  Therefore,  once  is  greater  than  or 
equal  to  Nq,  reliable  communication  is  drastically  reduced. 

For  the  switched  broadband  jamming  situation,  compare  the  switched 
broadband  jamming  equation  in  Figure  13  with  the  no  jamming  equation. 
Again,  they  are  of  the  same  basic  form.  However,  the  noise  in  the 
jamming  case  is  the  sum  of  and  K(m,n),  which  is  the  correlation 
coefficient  generated  because  the  duty  cycle  is  less  n  one  (see 

Equation  3-7).  Using  a  duty  cycle,  d,  equal  to  0.5, 


K(m,n) 


d=% 


J  /8 
o 


m=n 


m^n 
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Now,  the  noise  becomes  N  .  As  J  increases,  —  of  Figure  14 

ooo  N  b 

decreases  and  the  P^Ce)  increases.  The  same  result  as  in  the  broadband 

jamming  case  occurs.  Only  here  more  jamming  power  is  required  to 

E 

accomplish  the  identical  result.  That  is,  let  10  Log^  equal  13  dB 

o 

as  before.  In  order  for  this  ratio  to  drop  3  dB,  Jq  must  equal  4Nq. 
Once  Jq  surpasses  4Nq,  reliable  communication  is  reduced.  The 
modulator  is  operating  below  FM  threshold. 

Turning  to  the  continuous  wave,  or  single-tone  jamming,  the  CW 
jamming  equation  of  Figure  13  is  compared  to  the  no  jamming  equation. 
These  P^(e)  equations  differ  in  one  term  where  the  jamming  signal 
corresponds  to  one  of  the  eight  frequencies  of  the  modulator.  In  the 


I 

\ 

I 

n 
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demodulator,  the  received  signal  coefficients  are  pairwise  compared 
to  find  the  largest  one.  Without  jamming,  Equation  2-5  shows  that 
one  of  the  eight  coefficients  has  a  greater  expected  value  (mean)  than 
the  other  seven.  With  this  CW  jamming  the  mean  of  one  of  these  smaller 
coefficients  increases  such  that  the  probability  distribution  function 
of  Equation  2-5  becomes 


fr 

m 1  l 


N(/F+  j  ,  4N  ) 
m  o 

m=i 

N(  j  4N  ) 

m  o 

m^i  or  k 

4-2 

n(VTT,  sn  ) 

0  0 

m=k 

where  k  identifies  the  function  $>  (t)  which  has  the  same  frequency  as 
the  jammer  and  j  depends  upon  whether  the  jamming  frequency  equals  one 
of  the  signaling  frequencies  or  not.  Figure  15  illustrates  these 
probability  distribution  functions.  By  increasing  the  jamming  signal 
energy,  J  ,  the  mean  of  r  increases.  This  in  turn  increases  the 

O  K 

probability  of  selecting  r^  over  rm  as  the  largest  coefficient.  Let 
r^  equal  A,  as  shown  in  Figure  15.  With  a  mean  proportional  to  /J  , 


it  is  easy  to  see  that 


P[rk>A]  >  P[rn>A] 

So  r^  has  a  better  chance  of  being  the  largest  coefficient  than  any  of 
the  other  six  incorrect  coefficients  and  thus  disrupts  the  decision 
process  by  forcing  the  demodulator  to  return  the  wrong  channel  symbol. 

The  term  in  the  CW  jamming  equation  of  Figure  13  that  identifies  with 
r^  is  the  term  P[rk<A].  Remember  P^O)  =  — [1-P(C)J.  The  probability 
of  a  correct  decision  relies  on  the  product  of  the  pairwise  probabilities. 
The  pairwise  probability  for  r^  is  P [ r ^ <A ] .  So  as  J  increases,  from 
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the  jammer's  point  of  view,  P[r^<A]  decreases  which  in  turn  increases 
the  .  If  Jq  were  to  increase  such  that  h  / JqT  >  /E~,  then  the 

expected  demodulator  output  would  be  r^  most  of  the  time. 

Multitone  jamming  is  very  similar  to  CW  jamming  except  it  has 
more  than  one  jamming  tone.  In  the  equation  comparison  from  Figure  13, 
the  difference  between  the  no  jamming  P^Ce)  and  the  multitone  P^Ce) 
are  the  k  terms  which  identify  with  the  following  equation: 


k 

Tr  Pfr.<A] 
i-1  1 


4-3 


where  these  r^'s  correspond  to  4>^(t)'s,  which  in  turn  have  the  same 
frequencies  as  the  jammer.  Again  as  Jq  increases  this  product 
decreases.  As  a  result  P^Ce)  increases;  in  this  case  it  would  increase 

E 

more  rapidly  than  with  one  tone.  Should  J  become  greater  than  4  s, 

o 

then  the  k  jamming  signals  would  dominate  the  decision  process  of  the 
modulator.  The  same  three  probability  functions  discussed  previously 
with  CW  jamming  apply  for  the  multitone  jamming  case.  The  difference 
lies  in  the  number  of  coefficients  described  by  each  distribution.  With 
k  tones,  f  (x)  applies  to  each  one.  This  makes  more  coefficients 
more  likely  to  become  greater  than  r  . 


Encoder-Decoder  Bit  Error  Probabilities 

Appendix  B  contains  the  mathematical  process  of  bounding  the  bit 
error  probabilities.  These  two  probabilities  are  different  from  those 
in  Figure  13.  The  bit  error  probability,  P^,  is  defined  as  the  expected 
number  of  bit  errors  in  a  given  sequence  of  received  bits  normalized  by 
the  total  number  of  bits  in  the  sequence  for  particular  coders.  The 
results  of  the  bounding  for  each  coder  is: 
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CC1 : 


Dual-3: 


where  T(D,I)  is  a  generating  function  which  identifies  the  different 

Hamming  distances  between  an  incorrect  and  correct  codeword  and  the 

number  of  decoder  output  bits  that  are  in  error  (28:248).  Appendix  B 

identifies  T(D,I)  for  the  dual-3  decoder  error  events  which  are  two 

or  three  branches  long  in  the  trellis  of  Figure  B-4.  The  determination 

E 

of  larger  error  events  would  not  aid  in  the  analysis,  since  D  =  exp  -  rp-. 

Es  . 

As  —  increases,  D  approaches  0,  and  as  this  ratio  decreases  D  approaches 
o 

1.  Figure  16  shows  this  exponential  relationship  between  D  and  this 

E 

signal-to-noise  ratio.  For  ^  =  10,  D  is  essentially  0,  and  the  Pb(e) 

‘  o 

is  also  0  for  both  cases. 

Figure  17  shows  the  results  of  bit  error  rate  simulations  conducted 

by  Viterbi  (9:840)  for  a  soft,  8-level,  receiver  quantization.  Even 

though  he  restricted  his  simulation  to  32-bit  paths,  his  results  closely 

approximated  the  upper  bounds.  (For  low  error  probabilities,  i.e., 

Eb 

—  >  5  dB,  the  upper  bound  lies  slightly  below  the  simulation.  However, 
o 

the  quantization  loss  is  on  the  order  of  0.25  dB  which  appears  to  be  the 
approximate  separation  between  his  simulation  and  the  upper  bounds  (28:248).) 

With  the  above  data,  consider  the  broadband  jamming  case,  and  its 
effect  on  the  signal-to-noise  ratio.  The  system  noise  is  now  the  sum  of 
the  channel  AWGN  Nq  and  the  jamming  AWGN  J  .  Referring  to  the  upper 
bound  equations  of  P^,  as  Jq  increases  the  Q(-)  function  and  the  exponen¬ 
tial  function  of  T(D,I)  increase.  The  graph  of  Figure  17  better  depicts 
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CC1  SIMULATION  DATA 


FIGURE  17: 


P.  vs 


\r 


what  an  increase  in  channel  noise  does  to  the  communication  system.  As 

signal-to-noise  decreases  (moves  to  the  left),  the  P,  increases.  This 

b 

effect  also  holds  true  for  the  switched  broadband  jamming  case  except 
that,  with  a  duty  cycle  of  0.5,  four  times  as  much  jamming  energy  is 
required  to  accomplish  the  identical  effects  of  straight  broadband 
j  amming . 

The  effects  of  CW  jamming  are  probable  in  two  areas.  These  are 
noise  energy  per  transmitted  symbol  and  the  effects  of  phase  variations. 

For  the  CIV  jammer,  the  average  noise  power  affecting  the  receiver  is 
approximately  or  the  average  energy,  considered  as  jamming  noise, 
is  approximately  •  Letting  the  average  energy  of  the  jamming  signal 
be  J^, the  effective  noise  energy  of  the  channel  is  the  sum  of  and  . 

However,  this  now  becomes  nonwhite  noise  interference.  In  the  presence 
of  nonwhite  noise  an  optimum  receiver  can  be  designed.  However,  the 
effects  on  an  optimum  Gaussian  noise  receiver  are  in  the  spectrum  of 
the  received  signal.  The  nonlinearities  would  increase  the  noise  in 
the  system,  and  the  basic  effect  would  be  to  decrease  the  signal-to-noise 
ratio  and  thus  increase  the  probability  of  bit  error. 

The  second  area  of  probable  cause  for  system  degradation  would  be 
in  the  area  of  phase.  The  jamming  signal  has  a  phase  component  (see 
Equation  3-11).  In  the  coherent  detection  receiver  know  the  phase  of 
a  carrier  in  an  FM  system.  Receiver  component  instabilities,  uncompensated 
signal  shifts  and  noisy  signals  will  allow  less  than  optimum  phase 
reference  in  some  receivers.  In  a  correctly  received  signal,  its  energy 
component  can  be  reduced  by  a  factor  of  Cos  6  where  0  is  the  phase  shift 
(9:844).  Now  suppose  the  CW  jamming  signal  matches  the  one  of  the  eight 
signals  that  is  carrying  the  information.  With  the  proper  phase 
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manipulation  this  received  signal -to-noise  ratio  (S/M)  could  be  reduced, 
i  .e .  : 


c  E  J 

S  S  G 

si =•  «r  4  j-  cos  9 
'  o  '  o 


4-4 


As  -  approaches  this  S/M  approaches  0.  The  effect  of  phase  problems 
on  P^  is  critical  because  of  the  steep  bit  error  probability  versus 
signal-to-noise  ratio  curves  in  Figure  17.  A  slight  loss  of  may 
increase  drastically. 

Multitone  jamming  is  an  extension  of  the  CW  jamming  problem.  The 
effects  of  multitone  jamming  could  only  amplify  the  effects  of  nonwhite 
noise  and  phase  variations.  As  more  of  the  bandwidth  of  the  receiver 
is.  filled  with  multitone  jamming  signal  frequencies,  the  effective 
nonwhite  noise  approaches  white  noise  with  infinite  jamming  frequencies. 
With  each  of  the  tones,  an  average  energy  value  can  be  obtained.  The 
sum  of  all  these  average  energy  values  then  present  the  nonwhite  noise 
energy  of  the  system.  Again  its  effect  on  the  communication  system  is 
proportional  to  the  use  of  the  receiver  spectrum.  These  nonlinearities 
would  increase  the  noise  component  in  the  signal-to-noise  ratio.  The 
graph  in  Figure  17  illustrates  that  P^  increases  when  the  noise  increases. 

The  second  area  of  phase  variations  increases  the  effects  from  the 
CW  jamming  case  since  there  are  more  jamming  tones  available  to  enhance 
and  degrade  received  signal  energy.  Let  J  be  the  average  energy 
available  from  the  jammer  at  a  frequency  equal  to  the  frequency  of  S.(t). 


With 


a  proper  setting  of  the  phase  angle,  0,  S^(t)  can  be  made  to  appear 


as  the  transmitted  signal  if  Je^=Es-  In  the  same  light,  if  S^(t)  were 
the  transmitted  signal,  then  any  one  of  the  jamming  s  gnuls  with  the 
right  phase  and  energy  could  cancel  the  signal's  effect  in  the  receiver. 


56 


Conclusions  in  Performance  Analv sis 

This  Chapter  addressed  the  effects  of  four  types  of  jamming  on 
the  channel  model  characterized  in  Chapter  II.  This  performance  analysis 
covered  the  effects  on  the  probability  of  bit  error  and  the  coding  scheme 
bit  error  probability. 

The  broadband  jammer  affects  the  signal -to-noise  ratio  directly. 

As  an  AWGN  source,  the  noise  of  the  system  is  the  sum  of  .J  and  N'  . 

o  o 

With  increased  jamming  signal  energy,  the  signal - to-noise  ratio 

decreases.  The  corresponding  increase  in  probability  of  bit  error,  P 

can  be  followed  with  the  change  of  s ignal -to-noise  ratio  in  Figure  14. 

The  expected  values  of  would  coincide  with  the  graph. 

In  the  coding  bit  error  rates,  as  the  signal -to-noise  ratio 

decreases  the  bound  on  the  actual  probabilities  increases  for  both 

coders.  Figure  17  shows  this  relationship. 

Switched  broadband  jamming  has  the  same  effect  on  the  channel 

model  signal -to-noise  ratio  as  the  broadband  c  .-,e ,  except  with  a  duty 

cycle  of  0.5  the  noise  component  is  only  '..J  .  However,  it  increases 

P^(e)  in  accordance  with  Figure  14.  It  also  increases  the  bit  error 

probability  bounds  as  shown  in  Figure  1T. 

Continuous  wave  jamming  presents  a  different  effect  on  the  channel 

model.  It  principally  inhibits  the  demodulator's  ability  to  make  a 

correct  decision.  It  affects  the  pairwise  comparison  of  the  received 

coefficients  by  increasing  the  probability  that  one  of  the  seven  wrong 

coefficients  will  be  greater  than  the  correct  coefficient.  Figure  15 

shows  one  of  the  coefficients  having  a  different  mean  based  upon  the 

4F. 

s 

energy  in  the  jamming  signal.  If  J  >  — then  the  coefficient  associated 


with  that  jamming  signal  will  become  the  expected  output  of  the  demodulator 


The  effect  CIV  jamming  has  on  the  coders'  bit  error  rate  is  twofold. 
First,  by  virtue  of  the  signal  not  being  part  of  the  message  signal, 
it  adds  noise  to  the  channel.  This  additional  noise  decreases  the 
signal-to-noise  ratio  and,  thus,  increases  the  bit  error  probability. 
Second,  the  phase  of  the  jamming  signal  may  cancel  the  intended 
transmitter  signal  or  enhance  another  signal  to  pose  a  coding  problem 
of  two  signals.  In  either  event  the  signal-to-noise  ratio  is  decreased, 
and  Figure  17  shows  this  relates  to  an  increase  in  the  bit  error 
probability. 

Multitone  jamming  is  similar  to  CW  jamming  except  there  are  more 

tones.  An  increase  in  the  probability  of  the  pairwise  comparisons  being 

wrong  results  in  an  increase  in  the  P^(e)  .  Figure  15  shows  the  result 

of  increasing  the  expected  value  of  the  incorrect  coefficients.  If 

4E 

the  jamming  energy,  J  for  each  tone  surpasses  -j-  then  those  are  the 
coefficients  that  will  dominate  the  demodulator  output.  If  the  correct 
coefficient  is  not  one  of  these  coefficients  then  this  jamming  has 
successfully  blocked  the  communication  system. 

Just  like  the  CW  jamming  case,  multitone  jamming  can  affect  the 
bit  error  probabilities  of  the  coders  by  presenting  noise  by  their 
existence  to  the  channel  model.  This  noise,  being  nonwhite,  will 
decrease  the  signal-to-noise  ratio  and  increase  the  bit  error 
probability.  It  can  also  present  enough  phase  variations  to 
significantly  enhance  signals  not  present  and  degrade  the  transmitted 
signal . 
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V.  SIMULATION  PROGRAM  AND  SUBROUTINES 


Introduction 

This  Chapter  describes  the  simulation  program  for  the  communication 
system  in  Figure  9.  The  following  discussion  will  identify  the  basic 
program  elements,  describe  how  they  interact,  and  discuss  program 
validation,  problems,  and  results.  This  program  is  coded  in  FORTRAN  V 
for  the  CDC  6600-class  machine  belonging  to  the  Avionics  Systems 
Division  of  Air  Force  Systems  Command.  FORTRAN  V  was  selected  for  two 
reasons:  1)  effective  string  manipulations  are  possible  and  2)  the 

simulation  of  this  communication  system  is  not  time  dependent.  Such 
languages  as  GPSS  or  Simula  require  a  time  dimension  in  their  simulations. 
A  complete  program  listing  is  in  Appendix  C. 

The  organization  of  this  Chapter  begins  with  a  few  assumptions 
needed  for  programming  simplicity.  These  assumptions  create  no  loss 
in  generality  nor  deviation  from  the  thesis  objective.  Next,  Figure  IS 
shows  how  the  main  program  is  divided  into  its  major  operating  blocks. 

The  initialization  paragraph  discusses  parameter  declarations  and 
initial  parameter  values.  The  Encoder  section  identifies  the  two 
encoders  and  describes  the  overall  coding  operation  as  bits  in  and 
codewords  out.  The  next  part  of  the  main  program  is  the  Modulator. 

In  this  paragraph  the  modulation  scheme  is  described  and  explained  in 
terms  of  program  variables.  The  Jammer  section  selects  and  calls  the 
particular  jamming  subroutines.  The  last  section  of  the  main  program 
discusses  the  operation  of  the  Demodulator  and  the  call  of  the  decoding 
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subroutines . 


The  description  of  the  subroutines  falls  into  three  categories: 
the  Encoder,  the  Jammer,  and  the  Decoder  subroutines.  Each  subroutine 
has  a  flowchart,  identifies  its  function,  and  discusses  its  implementation 
The  eight  subroutines  are: 


ENCDR1 

ENCDR2 

BBJM 

SBJM 

CWJM 

MWJM 

VITDC1 

VITDC2 


Encoder  Subroutines 


Jammer  Subroutines 


Decoder  Subroutines 


The  final  portion  of  this  Chapter  discusses  the  program  validation 
and  simulation  results.  Expected  error  probabilities  are  compared  with 
empirical  values. 

In  order  to  build  this  simulation  model  certain  assumptions  were 
necessary.  These  assumptions  are  the  first  section  of  this  Chapter. 


Assumptions 

These  basic  assumptions  apply  to  the  coding  and  operation  of  this 
program. 

1.  The  input  data  can  be  either  a  "1"  or  a  ”0".  However,  since 
each  bit  is  equally  likely,  the  input  array,  INDATA,  is  initialized 
to  contain  all  0's.  With  this  input,  no  loss  in  generality  occurs,  and 
a  "1"  at  the  decoder  output  constitutes  an  error. 


2.  The  input  block  is  limited  to  100  bits. 

3.  The  signal  energy  is  unity.  This  preserves  the  orthonormality 


of  the  signals  and  simplifies  the  calculations.  Thus,  the  signal-to-noise 
ratio  is  the  inverse  of  the  noise  power  parameter,  NO. 


4.  The  encoders  are  assumed  to  start  and  stop  in  the  0  state,  i.e., 
00  state  for  encoder  1  and  000  state  of  encoder  2. 

5.  Synchronization  is  assumed  at  all  times.  This  means  that  the 
beginning  and  ending  bits  of  the  codewords  are  known.  No  markers  or 
pulses  are  required  to  obtain  synchronization.  Assuming  synchronization 
were  a  problem,  the  correct  path  through  the  trellis  will  not  emerge  as 
the  path  with  the  smallest  Hamming  distance.  In  fact,  no  path  will 
emerge  as  the  best  one.  All  cumulative  Hamming  distances  will  remain 
relatively  close  together.  Once  this  condition  is  identified,  the 
original  starting  bit  can  be  changed  to  an  adjacent  bit.  The  number 

of  starting  bit  changes,  or  shifts,  will  be  no  greater  than  the  codeword 
length,  n.  Synchronization  can  then  be  achieved  within  a  few  hundred 
bits  when  n  is  on  the  order  of  four  (28:261). 

Main  Program 

The  main  program  is  organized  into  functional  areas  of  the 
communication  system  in  Figure  9.  These  areas  are  the  subjects  of 
the  following  sections.  The  first  of  these  sections  is  Initialization. 
Initialization .  To  begin  the  program,  all  variables  and  arrays  are 
declared  according  to  type  and  initialized.  The  program  comments 
identify  each  variable  and  array.  The  initialization  of  parameters 
follows  the  variable  declarations.  Beginning  at  program  statement 
label  50,  the  random  number  subroutines  have  their  seeds  initialized. 

Two  library  procedures  are  used  to  generate  random  numbers.  The 
variable  SD4  is  the  seed  for  the  IMSL  uniform  random  generator,  GGUBFS(-). 
The  other  three  variables:  SD1,  SD2,  and  SD3,  are  used  in  the  normal 


random  number  generator,  GGNQF(-),  in  the  AKGN  calculation  and  two 
jamming  subroutines.  The  initial  values  of  these  three  seeds,  when 
used  in  parallel  with  the  single  random  generator,  create  uncorrelated 
random  number  strings.  The  seeds  are  initialized  by  the  same  number 
to  allow  for  comparisons  of  output  data  between  different  input 
parameters . 

Parameter  and  program  variations  are  controlled  by  use  of  the 
READ  statement.  Comments  describe  each  input  variable.  The  READ 
statement  has  a  free  field  format  for  simplicity.  The  Echocheck  is 
immediately  after  the  READ  statement  and  verifies  that  the  values  stored 
in  the  input  variables  are  the  same  as  the  ones  read  from  the  data  cards. 
The  last  portion  of  the  initialization  process  is  the  recognition  that 
the  source  output  bit  stream  (...  b^,  b^  ,  b^  ?,  ...)  contains  all 
zeroes.  Since  the  probability  of  a  ”1"  or  a  ”0"  is  the  same,  by  using 
all  zeroes,  a  "1"  at  the  destination  is  an  error. 

Encoders .  The  encoder  block  of  this  program  takes  the  string  of 
zeroes  and  sends  them  to  one  of  two  encoder  subroutines.  The  choice  of 
encoder  is  based  upon  the  value  of  the  input  variable,  ICC.  If  ICC 
equals  "1"  then  the  CC1  encoder  (ENCDR1)  is  called;  otherwise,  the 
dual-three  encoder  (ENCDR2)  is  called.  The  only  restriction  in  calling 
these  routines  is  that  100  bits  are  used  per  input  data  block.  These 
subroutines  are  explained  later  in  the  Encoder  Subroutines  section. 

From  the  encoders,  the  codewords  are  stored  in  the  C  array  and  are 
passed  to  the  Modulator. 

Modulator.  Figure  19  shows  the  8-ary  FSK  modulator  and  how  the  signals 
are  generated.  The  C  array  elements  are  grouped  into  actual  digits 
(0,  1,  ...  7)  and  stored  in  the  transmission  signal  array,  S(-).  The 


variable  K,  which  is  equal  to  70,  is  the  total  number  of  transmission 
signals  created  from  the  C  array.  Each  octal  digit  represents  one  of 
the  actual  modulator  message  signals  discussed  in  Chapter  II.  (See 
Equation  2-2.)  For  this  program  only  the  octal  digit  is  necessary  to 
simulate  the  transmission  signal  and  retain  the  message  information, 
.hammer.  In  this  section  of  the  main  program,  the  variable  JAMMER 
specifies  which  jamming  subroutine  will  be  used.  JAMMER  is  set  in  the 
READ  statement.  For  a  value  of  1,  the  broadband  jamming  subroutine, 

BBJM,  is  called;  for  a  value  of  2,  the  switched  broadband  jamming 
subroutine  is  called.  Likewise,  for  values  of  3  or  4,  the  continuous 
wave  jamming  (CWJM)  and  the  multi  tone  jamming  (MWJM-'  subroutines  are 
called,  respectively.  Each  jamming  subroutine  is  discussed  later  in 
the  Jammer  Subroutines  paragraphs.  A  series  of  IF -THEN -ELSE IF 
statements  insure  that  only  one  jamming  subroutine  is  selected.  If 
JAMMER  does  not  equal  one  of  these  four  numbers,  then  the  program 
assumes  no  jamming  is  present  and  proceeds  with  the  simulation. 
Demodulator.  This  portion  of  the  program  calculates  the  R(J)  co¬ 
efficients  and  determines  which  one  is  the  largest.  Within  the  trans¬ 
mission  medium  the  transmission  signal  is  affected  by  the  AWGN  of  the 
channel  and  the  jamming  signal,  if  present.  To  simulate  this  interference 
effect,  each  received  coefficient  is  calculated  based  upon  the  transmitted 
message  signal,  the  AWGN,  and  the  jamming  signal.  Figure  20  shows  the 
computation  of  the  received  coefficients  R(J) .  This  computation  is  the 
one  in  Equation  3-2,  allowing  the  array  JN(J)  to  represent  the  effects 
of  any  one  of  the  jamming  techniques.  The  integer  I  identifies  the 
transmitted  signal,  S(I);  NO  is  the  noise  power  parameter;  WN  is  the 
I  AWGN  coefficient;  GGNQF(-)  is  a  normally  distributed  random  number 
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FIGURE  19:  8- ARY  FSK  MODULATOR 

where  K  is  the  number  of  transmitted  signals 
and  M  is  the  Encoder  bit  memory  size. 


FIGURE  20:  RJ  COEFFICIENT  CALCULATIONS 


generator,  N(0,1);  JN  ( J )  is  the  jamming  signal  coefficient  array;  and 
RHO(-,-)  is  the  correlation  coefficient  matrix  of  a  set  of  orthonormal 
functions.  The  transmitted  signal  energy  is  1  to  preserve  the  ortho¬ 
normality  and  to  make  the  calculations  simple.  Each  time  a  coefficient, 
K(J),  is  calculated,  a  comparison  is  made  to  find  the  largest  R(J) . 

The  final  value  of  INDEX  is  the  subscript  of  the  largest  coefficient. 

The  demodulator  programming  block  transforms  INDEX  into  a  binary 
number,  three  bits  long.  The  range  of  values  for  INDEX  is  zero  to 
seven.  With  each  received  signal  coefficient,  a  three-bit  word  is 
concatenated  with  all  the  previous  three-bit  words  to  build  the  CH 
array.  This  array  is  the  estimate  of  the  C  array  which  contains  the 
coded  bits  entering  the  modulator.  Once  the  CH  array  is  filled,  it  is 
ready  for  decoding. 

Decoder.  One  of  two  subroutines,  each  using  the  Viterbi  algorithm 
(10:211),  processes  the  CH  array  and  constructs  the  output  bit  string 
in  B(0).  Array  element  B(0)  represents  the  state  in  which  the  encoder 
finished  its  coding  process  with  the  input  data  block.  More  detailed 
explanations  of  the  subroutine  operations  are  in  the  next  few  sections. 
The  final  operation  of  the  decoder  is  to  print  out  all  100  bits  obtained 
in  the  decoder  subroutines. 

Finally,  by  means  of  the  input  variable,  ISZ,  the  basic  program 
determines  whether  or  not  there  are  more  input  data  blocks.  If  so, 
the  program  returns  to  the  encoder  input  stage  and  begins  processing 
the  next  block  of  100  bits.  If  there  are  no  more  data  blocks,  the 
program  jumps  back  to  the  READ  statement  to  reinitialize  some  input 
variables  and  begin  processing  from  this  point. 

The  program  ends  upon  a  no-input-data  determination.  The  default 
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j ump  is  located  in  the  READ  statement  by  END=600.  When  the  input  data 
are  exhausted,  the  program  jumps  to  the  program  statement  label  600 
and  terminates  the  simulation. 


Desert ption  of  Subroutine s 

ENCDR1 .  This  encoding  subroutine  is  based  upon  the  description  of  the 
CC1  encoder  in  McEliece,  Chapter  9,  (11:200).  The  input  parameters  are 
the  block  size,  I,  and  the  input  data  .array,  IA.  There  are  three  outputs: 
the  coded  bit  stream,  C;  the  convolut ional  coding  transition  matrix,  CCM; 
and  the  number  of  bits  in  the  memory  register,  M.  Figure  21  shows  both 
the  matrix  and  the  coding  block  diagram.  Referring  back  to  Figure  2, 
the  explanation  of  the  transition  matrix  is  straight  forward.  Either  a 
”1"  or  a  ”0"  is  at  the  input,  thus  not  all  transitions  are  possible.  For 
instance,  given  a  0  present  state,  in  two-bit  representation:  00,  and  a 
"1"  input,  the  next  state  is  10  or  a  2.  Therefore,  a  1  state  (01)  will 
never  follow  a  0  state  (00).  The  output  bits  are  within  the  transition 
table.  For  continuity  within  the  program,  these  output  symbols  are  octal 
numbers.  Thus,  each  element  of  CCM  is  constructed  in  the  following  manner 


CCM(6)  = 


( _ i - - output  symbol 

3  1  0  0  1  0  0  0 

second  octal  number 
first  octal  number 

state 


This  element  represents  the  transition  from  the  third  state  (11)  to  the 
first  state  (01)  with  the  output  symbols:  001000.,  or  10g. 

Figure  2 a  also  shows  the  modulo-2  addition  which  calculates  the 


output  symbols:  C.,^  ^C,^.  Figure  21b  contains  the  subroutine  block 
diagram  for  this  encoder. 


ENCDR2 .  This  encoder  came  from  the  Mission  Research  Corporation  report 
(2:42)  and  is  known  as  a  dual-3  encoder.  The  inputs  to  this  subroutine 
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are  the  same  as  for  EXCDR1 :  I  is  the  input  block  size  and  IA  is  the 
input  data  array.  The  outputs  are  likewise  the  same:  C  is  the  coded 
bit  stream,  CCM  is  the  transition  matrix,  and  M  is  the  number  of  memory 
bits  in  the  encoder.  The  transition  table  is  already  shown  in  Figure  3b. 
Each  element  of  CCM  is  constructed  in  the  same  manner  as  for  EXCDR1 
above.  An  example  of  an  element  of  CCM  is  the  following  with  an 
explanation : 

CCM (35)  =43111010 

This  element  represents  the  transition  from  the  fourth  state  (100)  to 
the  third  state  (Oil)  with  the  output  symbols  111010.,  or  720.  All 
transitions  are  possible  by  the  nature  of  this  encoding  scheme. 

Figure  5a  shows  how  the  individual  bits  are  combined  to  form  a  coded 
symbol  and  how  they  are  shifted  for  each  new  three-bit  input.  Thus, 
all  possible  transitions  are  allowed.  The  flowchart  for  this  subroutine 
is  shown  in  Figure  22. 

Jammer  Subroutines 

BB.JM .  This  subroutine  represents  the  baseband  jamming  process.  The 
inputs  are  the  random  number  generator  seed,  SD2,  and  the  power  parameter 
in  the  jamming  signal,  JO.  The  output  is  the  array,  JN ,  which  contains 
the  eight  jamming  coefficients,  j  ,  from  Equation  3-2.  To  arrive  at 
each  jamming  coefficient,  JN(I),  with  a  distribution  of  N(0,'jJ0),  a 
random  number  is  obtained  from  the  IMSL  library  routine,  GGNQF(-)>  and 
multiplied  by  the  variance  of  JN'(I).  Eight  of  these  coefficients  are 
calculated  and  returned  to  the  main  program  to  be  used  in  the  eight 
receive  coefficients  in  the  demodulator. 

SB.JM .  This  subroutine,  for  the  switched  broadband  jammer,  is  more 
complicated  than  the  broadband  case.  Through  the  calculations  in 


Chapter  [II  and  Appendix  A,  the  eight  coefficients  are  normally 
distributed  with  the  variance  and  mean  identified  in  liquation  3-8. 

Figure  do  is  the  flowchart  for  this  subroutine.  The  three  inputs  are 
the  random  number  generator  seed,  SD5;  the  power  parameter  in  the 
jamming  signal,  .10;  am.  *he  average  duty  cycle,  DC.  The  output  is  the 
8-arv  array,  J\' ,  containing  the  eight  jamming  coefficients.  The 
variable  ALPHA  is  the  signal's  angular  frequency  used  to  calculate  the 
variance  of  the  jamming  random  variable.  Again  GG.\’QF(.)  is  used  with  a 
different  seed  from  before  to  finally  calculate  these  coefficients  which 


are  transferred  back  to  the  main  program. 

CWJM .  The  continuous  wave  jammer  subroutine  is  Howcharted  in  Figure  34 
The  four  inputs  are  the  random  number  generator  seed,  SD4 ;  the  power 
parameter  of  the  jamming  signal,  JO;  the  jamming  frequency,  OMEGAJ  or 


co j ;  and  the  transmission  time  interval,  T.  The  output  is  the  eight 

element  array,  JN .  For  comparison  purposes,  uj  ,  the  signal  frequency, 

and  to.  are  kept  in  Hertz  until  the  subroutine  determines  whether  thev 
J 

are  linearly  related.  No  matter  what  the  relationship  is  between  ,.o-  and 


U)£,  JN(I)  is  calculated  according  to  Equation  3-12,  where  the  phase  is  0 
This  subroutine  handles  one  jamming  signal.  If  there  is  more  than  one 


tone,  the  subroutine  MWJM,  described  in  the  next  section,  controls  the 


effects  of  multiple  tones  on  the  received  signal. 


MWJM.  This  is  the  multitone  jamming  subroutine.  The  basic  equation  to 
calculate  ,JN ,  the  jamming  coefficient  array,  is  Equation  3-16.  The 
inputs  to  this  subroutine  are  the  random  number  generator  seed,  SD4 ; 


the  number  of  tones,  NUM;  the  power  parameter,  JO;  the  transmission 
interval,  T;  and  the  jamming  tone  frequencies,  WJ ( • ) .  The  output  is  the 
IN  array  which  has  the  eight  jamming  coefficients.  This  subroutine  call 
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the  CWJM  subroutine  for  each  tone  to  calculate  the  contributions  of 
these  tones  to  the  overall  array  JN . 


Decoder  Subroutines 

VITDC1 .  The  Viterbi  algorithm  is  the  basic  procedure  for  this  sub¬ 
routine.  Figure  25  is  the  subroutine  flowchart.  The  inputs  are  the 
array  with  the  estimated  coded  symbols  in  binary  form,  CM;  the  number 
of  elements  in  the  CM  array,  K;  and  the  transition  matrix,  CCM.  The 
lone  output  is  the  bit  string  in  the  array  element  B(0).  This  element 
estimates  the  original  input  bit  string  from  the  source. 

Internal  variables  and  arrays  are  as  follows: 

BP  "scratch-pad"  character  array 

HMD  the  Hamming  distance  array 

D  the  cumulative  distance  array 

IFM  present  transition  state 


next  transition  state 


J,  CNT,  $  ICH  integer  counting  variables 


1A  5  IB 


temporary  storage  locations 


The  HMD  array  is  organized  to  have  its  first  subscript  represent  the 
present  state  and  to  have  its  second  subscript  represent  the  next  state, 
i.e.,  HMD(0,2)  equals  the  Hamming  distance  between  the  encoder  output 
for  the  transition  from  state  0  to  state  2  and  the  received  codeword. 

The  D  array  cumulatively  adds  the  Hamming  distances  to  each  node  of  the 
coding  trellis,  see  Figure  7. 

The  intrinsic  function  ICHAR(arg)  returns  the  position  number  of 
the  alphanumeric  argument  in  an  intrinsic  ASCII  ordering  list.  This 
ASCII  list  has  positions  16  to  25  occupied  by  the  numerical  characters 


0  to  9.  Recalling  that  the  first  character  of  an  element  in  CCM  is  the 
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present  state  of  a  potential  transition  and  the  second  character  of  an 
element  in  CCM  is  the  next  state,  IFM  and  ITO  are  set  to  the  proper 
state  numbers  after  subtracting  16  from  the  ICHAR  value.  IFM  and  ITO 
indicate  which  transition  the  subroutine  will  consider  in  calculating 
the  Hamming  distances  for  entry  in  HMD.  The  received  codeword 
(CH(ICH) ,CH(ICH+1))  is  compared  with  the  remaining  characters  in  the 
element  of  CCM. 

The  final  portion  of  this  subroutine  uses  the  fact  that  only  eight 
transitions  are  possible  from  one  depth  to  another  in  the  coding  trellis. 
However,  only  four  transitions  can  be  carried  to  the  next  trellis  depth. 

Since  each  state  has  two  transitions  entering  it,  this  subroutine  selects 
the  transition  with  the  least  cumulative  Hamming  distance.  This  operation 
is  illustrated  in  Figure  7b.  The  program  calculates  the  potential 
cumulative  Hamming  distance  for  both  transitions,  i.e.,  IA  and  IB.  The 
smaller  value  becomes  the  next  state's  cumulative  Hamming  distance,  and 
the  appropriate  bit  character,  "0”  or  ”1”,  is  added  to  the  previous  bit 
string,  BP(-).  This  selection  process  is  done  for  each  state  and  each 
codeword . 

By  initializing  D(0,0)  to  0  and  the  other  first  row  elements  to 
999,  the  subroutine  automatically  assumes  that  the  starting  state  of 
the  encoder  was  0.  At  the  output,  the  first  100  bits  assembled  in  B(C) 
are  the  estimated  input  bit  string.  The  remaining  two  bits  are  a  result 
of  the  memory  bits  in  the  encoder.  Thus,  the  substring,  B(0)  (6:105), 
carries  the  output  bit  string. 

VITDC2 .  This  second  Viterbi  subroutine  uses  the  basic  algorithm  tailored 
to  the  second  encoder  discussed  in  Chapter  II.  Figure  26  is  the  subroutine's 
flowchart.  The  inputs,  identical  to  those  of  VITDC1,  are  the  array  with 
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the  estimated  coded  symbols  in  binary  form,  CH;  the  number  of  elements 
in  the  CH  array,  k;  and  the  transition  matrix  from  ENCDR2,  CCM.  Again, 
the  lone  output  is  the  assembled  bit  string  in  the  array  element  B(0). 
This  bit  string  estimates  the  original  input  string  to  the  encoder. 

The  internal  variables  and  arrays  are  as  follows: 

BP  "scratch-pad"  character  array 

D  the  cumulative  Hamming  distance  array 

HMD  the  Hamming  distance  array 

CNT,  J,  I,  5  ICH  integer  counting  variables 
CSV  present  state  storage  location 

DM  codeword  storage  location  for  comparisons 

IFM  present  transition  state 

ITO  next  transition  state 

Recall  that  HMD(1,3)  means  that  this  location  contains  the  Hamming 
distance  associated  with  the  transition  from  state  1  to  state  3.  The 
D  array  cumulatively  adds  the  Hamming  distances  at  each  node  of  the 
coding  trellis.  Also,  the  intrinsic  function  ICHAR(arg)  returns  the 
position  of  the  argument  in  an  intrinsic  ASCII  ordering  list.  The 
numerical  characters,  0  through  9,  occupy  positions  16  through  25  in 
this  ASCII  list.  Therefore,  by  subtracting  16,  this  list  position 
number  becomes  the  desired  numeric  value.  The  use  of  the  character 
variable  DM  distinguishes  this  subroutine  from  the  previous  one.  DM 
gathers  the  six  bits  that  make  up  the  transmitted  codeword.  By 
comparing  DM  with  each  transition  element  of  CCM,  the  Hamming  distance 
for  those  transitions  are  calculated. 

Once  the  HMD(*,-)  array  is  filled,  the  subroutine  proceeds  to 
determine  the  transition  to  each  stat^  at  the  next  trellis  depth.  This 
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means  that  for  each  state  there  are  eight  possible  transitions.  This 
subroutine  selects  the  transition  which  gives  the  smallest  cumulative 
Hamming  distance.  Each  transition  then  identifies  the  input  that  would 
make  it  happen  in  the  coding  scheme  discussed  in  Chapter  II.  In  the 
case  of  this  subroutine,  CSV  saves  this  input  and  is  translated  into  a 
binary  character  string,  which  is  concatenated  with  the  previous  input 
binary  character  strings.  Ultimately,  the  final  bit  string  is  in  the 
substring,  8(0)  (6:105),  which  is  passed  to  the  main  program.  The  extra 
character  array,  BP(-).  saves  the  previously  assembled  bit  strings  for 
each  decoding  state.  This  allows  the  formation  of  the  actual  array 
B(.)  without  altering  any  values  needed  in  completing  the  other  seven 
bit  strings. 

Validation  of  the  Simulation  Program 

Before  using  this  simulation  program  with  the  four  types  of  jamming, 
the  program  must  be  validated  as  a  legitimate  representation  of  an  AWGN 
digital  channel.  Three  separate  program  runs  were  made  to  check  this 
simulation  model.  Data  for  the  probability  of  symbol  errors,  probability 
of  bit  errors,  and  the  bit  error  rate  were  collected  and  compared  with 
the  expected  data  from  Chapter  IV  and  Figures  14  and  18. 

To  find  the  probability  of  symbol  error,  the  simulation  program  was 
modified  to  count  the  number  of  times  the  parameter  INDEX  was  not  equal 
to  zero.  By  varying  the  noise  parameter,  and  thus  the  signal-to-noise 
ratio,  the  symbol  errors  were  counted.  Figure  27  shows  the  expected  curve 
of  values  for  the  probability  of  symbol  errors  (29:259)  and  the  plotted 
simulation  points  (X)  for  various  signal-to-noise  ratios.  There  is  close 
agreement  between  the  expected  value  and  the  empirical  data. 

In  finding  the  probability  of  bit  error,  the  simulation  program  was 


78 


modified  this  time  to  count  the  non-zero  bits  placed  in  the  CH  array 
by  the  demodulator.  By  varying  the  signal-to-noise  ratio  the  bit 
errors  were  counted.  The  same  figure,  Figure  27,  shows  the  close 
agreement  between  the  expected  and  the  empirical  data  (+) . 

To  check  the  bit  error  rates  of  this  channel  mode,  data  was  collected 
from  simulation  runs  varying  the  signai-to-noise  ratio  for  each  encoder. 

The  entire  program  was  used  for  this  check;  that  is  to  say,  the 
simulation  was  from  the  encoder's  input  to  the  decoder's  output,  assuming 
AWGN  as  the  only  interference.  The  bit  error  rate  for  each  simulation 
was  the  ratio  of  the  number  of  ones  in  the  element  B(0)  to  the  total 
number  in  the  sample  size.  The  sample  sizes  were  only  1000  bits  because 
of  limited  computer  resources.  Figure  28  shows  the  comparison  for  each 
encoder  between  the  simulation  results  and  the  bounds  for  the  bit  error 
rates  found  in  Figure  18  of  Chapter  IV.  In  Figure  18  these  bounds  are 
over  a  wider  range  of  bit  error  rates.  Figure  28  illustrates  that  as 
the  bit  error  rate  decreases,  the  performance  of  the  dual-three  encoder 
is  better  than  the  performance  of  the  CC1  encoder. 

By  showing  that  this  simulation  program  models  and  AWGN  channel  in 
probability  of  errors  and  bit  error  rates,  the  model  is  now  ready  to  be 
tested  using  the  jamming  techniques  outlined  in  Chapter  III. 

Effects  of  the  Jamming 

Three  of  the  four  jamming  types  discussed  in  Chapter  III  were  applied 
to  this  simulation  program.  They  are  the  broadband  jamming,  switched 
broadband  jamming  and  the  continuous  wave. 

For  broadband  jamming  problem,  the  AWGN  was  assumed  to  be  part  of 
the  jamming  power  parameter,  JO.  Difficulties  arose  with  the  uncorre 1 ated- 
ness  of  the  same  random  number  generator  operating  in  parallel  with  itself 
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to  generate  n  and  j  .  Since  both  the  AWGN  nit)  and  the  jamming  signal 
J(t)  are  Gaussian  random  processes  in  the  channel,  their  sum  is  also  a 
Gaussian  random  variable.  By  letting  JO  contain  NO,  a  constant  for 
these  simulations,  the  signal -to-noise  ratio  varies  with  .It'.  For  this 
case  the  bit  error  rate  data  collected  for  channel  validation  were 
duplicated . 

For  the  switched  broadband  jamming  case,  the  same  results 
occurred  as  did  with  the  broadband  jamming  case.  Only  now,  ,J0  had  to 
be  four  times  as  great  to  produce  the  same  magnitude  of  change  as  it 
produced  in  the  broadband  case. 

The  CW  jamming  results  is  limited  to  testing  three  frequencies: 

900,  1800,  and  4800  Hertz.  Two  frequencies,  1800  and  4300,  are  from 
the  transmission  signal  set.  The  AIVGN  parameter  NO  was  0.10  for  this 
test.  Bit  error  rates  were  the  only  performance  measure  of  the  CW 
jamming.  Recalling  Figure  IS,  as  JO  increases,  the  mean  of  the  incorrect 
coefficient  increases.  Once  JO  is  greater  than  4Es/T,  the  bit  error 
rate  reflected  that  the  demodulator  output  of  this  signal  would  dominate 
the  bit  pattern.  For  instance,  using  the  CC1  encoder  and  un  equal  to 
1800  Hz,  the  expected  output  would  be  a  "011"  most  of  the  time,  when 
JO  >  4Es/T.  Decoding  a  string  of  "011"'s  gives  a  string  of  "100'"s  at 
the  decoder's  output.  As  JO  increased,  the  demodulation  and  decoding 
of  1800bps  would  dominate  the  output  and  the  bit  error  rate  would  tend 
towards  .333.  However,  if  the  same  type  of  reasoning  is  used  for  4800bps, 
the  the  bit  error  rate  would  tend  toward  1.  Therefore,  the  bit  error 
rates  are  dependent  upon  the  frequency  of  the  jamming  signal.  This 
dependence  is  illustrated  in  Figure  28.  The  data  for  each  jamming 
frequency  (i.e.,  1800  and  4800bps)  are  plotted  versus  the  expected 
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mean  of  the  incorrect  coefficient,  u.  =  !2(J  T)  .  As  u  increases  the 

J  0  j 

bit  error  rate  increases  as  expected.  The  reason  the  two  frequencies 
plot  different  curves  is  that  the  demodulator  and  decoder  translate 
each  one  into  a  different  bit  string. 
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V]  .  CONCLUSIONS  AND  RECOMMENDATION'S 


IHLscuss  ion 

Tiie  problem  stated  in  Chapter  I  identified  the  need  to  determine 
the  effects  that  ECM  has  on  error  correcting  _cdes  optimized  for  a 
non-ECM  environment.  The  narrowing  of  the  problem  to  two  convolutional 
encoders  and  four  jamming  models,  aided  in  building  a  specific  AWGN 
channel  with  specific  components.  However,  the  generality  of  the  problem 
statement  was  not  lost.  Now  that  the  channel  is  characterized,  any 
encoder-decoder  operation  or  jamming  technique  can  be  added  to  this  model. 
The  two  encoders  differed  greatly  from  simplicity  to  complexity  with 
even  their  coding  constraint  lengths  differing:  i.e.,  2  and  6.  The 
four  jamming  techniques  were  obtained  from  the  literature  search  and 
from  discussions  with  others.  Four  jammers  appeared  to  be  a  moderate 
number  to  use.  The  majority  of  work  for  this  thesis  centered  about 
two  aspects:  1)  the  research  and  understanding  of  convolutional  coding 
techniques,  Viterbi  decoding  algorithm,  and  the  development  of  the 
probability  of  error  terms;  2)  the  ’mplementat ion  of  the  simulation 
program  on  a  limited  computer  >y,--tem. 

Conclusions 

In  the  channel  performance  analysis  (Chapter  IV)  the  basic  conclusions 
are  based  upon  the  probability  of  bit  error  versus  signal -to-noi^e  ratio 
and  upon  the  bit  error  rates  for  each  of  the  encoder-decoder  operations.  In 
the  no  jamming  case,  it  appears  that  the  code  with  the  larger  constraint 
length  is  the  better  performing  code.  In  this  case  since  both  encoders 
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have  a  rate  of  4,  the  dual -three  encoder  performs  better  than  the  CC1 
encoder  (see  Figure  18).  With  an  additional  noise  source  present,  as 
indicated  in  Figure  9,  the  effects  on  the  coding  operation  varies. 

For  broadband  and  switched  broadband  jamming,  J(t)  is  a  Gaussian  random 
process  similar  to  the  AWGN  of  the  channel.  In  this  case,  both  encoders 
were  affected  similarly  and  the  dual-three  still  performed  better.  As 
with  the  two  previous  jammers,  the  continuous  wave  and  the  multitone 
jammers  affected  the  channel  model  similarly.  Depending  upon  where 
the  jamming  frequency  lies  in  the  transmission  bandwidth,  the  effects 
on  the  coding  performance  vary.  A  single  CW  tone  disrupts  the  decision 
process  of  the  demodulator,  which  in  turn  then  passes  an  incorrect 
message  codeword  to  the  decoder.  Since  the  dual -three  codeword  is 
longer  than  the  CC1  codeword,  the  decoding  process  of  the  dual-three 
codewords  is  the  more  susceptible  output  to  bit  errors  than  the  CC1 
decoding  process. 

The  simulation  program  (Appendix  C)  models  the  AWGN  channel.  By 
means  of  input  data,  either  encoder-decoder  operation  can  be  selected, 
as  well  as  one  of  the  four  jamming  subroutines.  This  program  was 
validated  by  checking  the  probability  of  bit  and  symbol  errors  and  the 
coding  bit  error  rate  against  known  data.  Chapter  V  discusses  the 
validation  process.  Simulation  of  broadband  and  switched  broadband  data 
resulted  in  the  same  data  as  for  the  channel  validation  where  the  AWGN 
parameter  was  varied  to  determine  the  system  response.  Figures  14  and 
18  apply  to  these  jamming  techniques.  The  CW  jamming  results  remain 
inconclusive  at  this  point  because  of  the  low  sample  size  collected. 
Appearances  indicate  that  the  bit  error  rate  is  affected  in  the  higher 
frequency  ranges  of  the  bandwidth  as  JO  becomes  greater  than  4Es/T. 
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This  is  expected  since  the  message  when  decoded  correctly  would  have 
no  l's.  Thus,  the  signals  representing  those  messages  with  l's  are 
higher  in  the  transmission  bandwidth.  Multitone  jamming  would  enhance 
this  last  point  by  giving  the  demodulator  a  choice  in  selecting  the 
wrong  message  signal  representation. 

Recommendations  and  Further  Study 

Since  computer  processing  time  became  a  premium  these  past  few 
months,  the  full  characterization  of  the  transmission  bandwidth  in  a 
CW  jamming  environment  was  not  possible.  Therefore,  the  first  step 
beyond  this  thesis  would  be  to  check  the  results  of  CW  jamming  when 
various  tones  from  the  bandwidth  are  used  for  the  CW  signal.  As  JO 
increases  the  bit  error  rate  should  stabilize  at  a  particular  value 
representative  of  the  decoding  of  the  jamming  signal. 

In  the  spirit  of  generality,  as  was  the  problem  statement  before 
the  problem  was  narrowed,  other  encoders  with  different  rates  and 
constraint  lengths  and  other  jamming  environments  can  be  made  to  fit 
in  this  simulation  model  and  analyzed  with  respect  to  one  another. 

Provided  better  computing  facilities  are  available,  the  further 
use  of  this  simulation  program  should  include  the  following  changes. 

The  actual  probability  of  symbol  and  bit  errors  and  the  bit  error  rates 
should  be  calculated  per  input  data  block.  This  will  save  manually 
counting  the  l's  and  0's.  Also,  the  number  of  samples  should  be  increased. 

_3 

With  1000  saemples,  the  best  bit  error  rate  is  10  .  Sample  sizes  should 

be  around  10  ^  to  be  able  to  compare  with  data  in  the  current  literature. 
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APPENDIX  A:  DERIVATION  OF  ERROR  PROBABILITIES 

This  Appendix  provides  the  background  equation  manipulations  and 
derivations  to  support  the  equations  in  Chapter  III,  Jamming  Models  and 
Performance  Analysis.  For  each  jamming  type  discussed  the  following 
information  is  basic: 


Input  Probabilities: 
P (b=l)  =  P(b=0) 

P(S  (t))  =  i 

1  a , 


=  h 


i 


1,  2, 


8 


Transmission  Signal  Set: 

Si(t)  =  / 2E/T  Sin(i2Trt/T) , 

te[0,T]  and  i  =  1,  2,  ...  8 

Modulator  Orthonormal  Signal  Set: 

^(0  =  (  |  )h  Sin(i2rr/T) ,  i  =  1,  ...  8 

AWGN  Probability  Distribution  Function: 

f  (x)  =  -i-  exp(- 

/  ttN  N 

o  o 

Jamming  process,  independent  of  the  AWGN  and  the  transmission 
signals,  is  J(t);  see  Figure  9  in  Chapter  III. 
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j  t 


Received  Signal: 


r(t)  =  Si(t)  +  n(t)  +  J(t),  te[0,T] 


Orthogonal  Coefficient  of  the  Received  Signal: 


rm  =/n  r(t)  WtJ  dt’  m  =  0>  1»  •••  7 


=  f  S  (t)  $  (t)dt  +  f  n(t)  $  .(t)dt  +  f  J(t)  $  .  (t)dt 

Jr \  1  m+1  Jr,  m+1  J  n  m+1 


=  /Fp.  +  n  +  f  J(t)  $  ,(t)dt 
m  Jq  m+1 


Broadband  Jamming 


T 

Let  j m  ~  f  J(t)  $m+i(t)dt.  ^en»  j's  Gaussian  Probability 
J  0 


Distribution  Function  is: 


f  (x)  =  — —  exp(~  J-)  a  N(0,5sJ  ) 

Jm  Jrrl  °  r> 


Therefore, 


r  =  /F  p.  ,  +  n  +  i 
m  i,m+l  m  i 


To  find  the  mean  of  r  ,  the  expected  value  of  r  is  found: 

m  m 


E{r  }  =  /F  p. 
m  Mi,m+1 


The  variance  of  r  is: 

m 


Var{r  }  =  E{r 2}  -  E2{r  } 
mm  m 


Therefore, 


r  2  =  Ep.  ,  2  +  n  2  +  i  2  +  fW  o ■  *(n  +  i  1 
m  Mi,m+12  m  Jm  Mi,m+1<'  m  Jnr 

•*  n  (/E~p.  ,  +  j  )  +  j  (/E~p.  ,  +  n  ) 

£{*  }  =  Epi  m+12  +  E^nm2^  +  E{j_2) 
m  i ,  m+ 1  m  J  m 


=  Ep.  2  +  ^  +  hJ 

^l.m+l*  o  o 


Var{r  }=  %(N  +  J  ) 

m  oo 


The  probability  distribution  function  of  r^  given  a  particular  hypothesis 
is  true  follows: 


fr  |Ht(x)  " 
m 1  k 


N(/F,  JjCN  +  J  )) 
o  o 

N  (  0,  *5(No  +  JQ)) 


m=k 

m^k 


When  S^+1(t)  is  transmitted  we  have: 


rk  =  /F  +  nk  +  jk 


r.  =  n.  +  i . 
i  i  i 


i*k 


Thus  the  probability  of  the  correct  r^  being  selected  as  the  largest  is 


P(C|Hk>rk)  -  P(r0<rk,  r1<rk.  ...  r7<rk) 


m^k 


The  last  equality  stems  from  the  fact  that  all  (n^  +  j ^) ? s  are  statistically 

independent  and  identically  distributed.  Multiplying  by  f  (x)  and 

rk  1  k 

integrating  gives: 


< 
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where 


Ptc|[v-  [v.:£rj «t ** 


P(r  <r,  ) 
m  k 


■/ 


r  _r  2 

k  1  r  m  ^  , 

—  exP  Cr; — — H  dr 
0°  /(No  +  JoJt;  No  +  Jo  m 


m^k 


let 


a  = 


m 


A(N0  *  Jo) 


da 


M  dr 

N  +  J  m 
o  o 


x  = 


v^Tn  +”  J  ) 
o  0 


dx  = 


N  +  J  drk 
o  o 


P(r  <r  ) 
m  k 


-i 


—  exp(-  J^a2)  da  =  l-Q(x) 
/2tt 


therefore. 


OO 

P(C|hk)  .  //xk  |Hk(V/r)  U-Q(x)]7  drk 


From  the  symmetry  of  the  model: 


7 

P(C)  =  Z  P(Hk)  P(C|H  )  =  P(C|H  ) 


dx 


00 

P(C)  =  -  f  exp 

/2tT 


P(e)  =  1  -  P(C) 

Switched  Broadband  Jamming 

The  jamming  process,  J(t),  is  the  same  as  in  the  broadband  case, 
except  there  is  now  a  switching  function,  z(t),  that  turns  this  jammer 
on  and  off.  The  following  information  characterizes  z(t)  with  an  average 
duty  cycle  of  d: 


P(z(t)  =  1)  =  d 

P(zCt)  =  0)  =  1-d 

E{z(t) }  =  d 

E{z2 (t) }  =  d 

Var  { z (t) }  =  d(l-d) 

J(t)  and  z(t)  are  independent;  therefore. 


r 

m 


Vi(t)  dt 


Let 


jm  ’  /„z(t)  J(t)  Vl(t)  dt 


■/ 
J  rv 


E(jm}  -  E{  /  2(t)  J(t)  Vl(.)  dt} 


E{z(t) }  E{j(t)J  $m+1(t)  dt 


0 


rT  rT 

E{-jmJk}  =  E{  J  Z(U^J(U)  $  ,(u)du  /  2(v)JCv)  $  (v)dv) 

'  r>  m  1  J  r\  k+l 


=  f  j1  E{z(u)z(v)}E{J(u)J(v)}  $  (u)  $,  (v)  dudv 

J  0  J  0  m+1  k+1 


T  T 

6(u-v)E{z (u) z (v) }  $m+1(u)  $k+1(v)  dudv 

rT 

o 


*  hJQ  f  E{z(u)z(u)}  $m+1Cu)  $k+1Cu)  du 


Since  z(t)  =  1  over  [0,dT]  and  0  elsewhere. 


dT 

E{jmjm}  =  H>d  /Q  Vl(u)  *k+l(u)  du 

Let  m=k  which  means 


E{jmjk}  =  E{V} 


dT  2 

EV,*'*Jod  /„  Vl  (t)  dt 


0 

dT 


=  ^  d  f  y  Sin2  ((m+l)47rt/T)  dt 
J  0 

=  ^od  [d  "  4rr(m+ 1)  Sin  ( (m+ 1)  4trdD  ] 


Now  let  m/k. 


dT 
( „ 
d.T 


E{jmjk> =  ^0d  /.  Vi(t)  dt 


=  5sJQd  J  |  Sin((m+l)4TTt/T)Sin((k+l)47Tt/T)dt 
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=  *od  1  27^10  Sin^2^m-k)d)-  2^i7k+2)  Sin(27r(m+k+2)d)] 


Let  K.(m,kj  be  the  covariance  function  for  i  . 
.1  m 


Vm>k]  =  E{jmV  *  E{V  E{V 


vod  [  d  -  4^iy  sin^m+1)«] 


H>d  i  sin(2T(n,'k)d) 


|  '  27(m+k+2)  Sin(2TT(m+k,2)d)] 

Since  J(t)  is  a  randon  Gaussian  process,  then  z(t)J(t)  is  a  random 
Gaussian  process.  The  probability  distribution  function  of  given 
then  becomes  the  following: 


“r  K 
m 1  k 


N(/E,  5sN  +  K .  (m,k) ) 
(x)  =  °  J 

N (  0,  lsNo  +  K.  (m,k)) 


Let  gZ=J5N0  +  K^(m,k).  As  with  the  broadband  jamming,  the  P(C)  is 
calculated  for  the  switched  broadband  jamming  similarly. 


p(C|Hk,rm)  -  P(ro<rk,r1<rk,  ...  r^) 


*  pcwi 


p(c|Hk}  =  //rk|Hk(x)  tP(rm<rk^  dx 
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From  the  model  symmetry  and  the  broadband  jamming  calculations: 
P(C)  =  P(C|Hk) 


l  r 30  2 

-  /  exp  [-  h(x  -  /I/jp!  ]  [—  e'^da]  dx 

/  2m  -00  /2  7T  -00 

where 

rm  rk 

a  =  -  ,  x  =  -  ,  g2  =  +  K.(m,k) 

_  _  o  j 

/g?  /g2 

P(e)  =  l-P(C) 


Continuous  Wave  Jamming 

Continuous  wave  (CW)  jamming  is  a  narrowband  process  where  J(t) 
represents  the  jamming  signal: 

J(t)  =  y/hJ0  Sin(w^t  +  0) 

Wj  is  the  jamming  frequency  and  8  is  the  phase  angle.  The  received 
signal  coefficierc  equals: 


r 

m 


$m+l(t)  dt 


T  /— 

=  ^Pi.m.!  +  nm  +/0  ^o  Sin(wjt  +  9)  A  Sin(wm°  dt 


where  wm  is  the  message  frequency  equal  to  (m+1) 2tt/T. 

Let  w.  ^  w  ,  then: 

J  m 


4>  ,(t)  dt 

m+  i 


Am 
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rv 

=  I  /j  /T  Sin(w.t  +  9)  Sin  (w  t)  dt 
J  n  °  3  m 


=  /J  /T  [  -  ][  4Cos9  Sin(aj.T)  +  4Sin0  (Costa. T)  -  1)] 


For  u).  ~  u  , 
J  m 


j  =  f  /J  /T  Sin  (w.t  +  0)  SinCto.t)  dt 

m  J  n  0  J  3 


/J  /T  4Cos0  T  =  /Tf  4Cos0 

O  0 


The  phase  angle  is  uniformly  distributed  over  For  this  paper, 
since  throughout  the  problem  coherent  detection  has  been  assumed,  the 
phase  is  0,  its  mean  value.  Therefore, 


*  *v 

2u> 

*5  SJir  [  — ■ T  J  Sin(u.T) 

°  J 


3  m 


3  m 


The  probability  distribution  function  for  rm  given  is: 


N(  /E  +  j  ,  %N  ) 

r  /  \  mo 

^r  Ih  = 

rm'Hk  N (  jmANQ) 


where  depends  upon  whether  w  equals  oj^  or  not.  The  probability  of 
correctness  and  error  follows  from  the  model  symmetry  and  the  previous 
calculations: 


P(C)  =  P(C[Hk) 


r  1 

/  T=exp 

An 

O 


frk  - (/f+  v)! 


[rANo  e> 
J-oo  /2tt 


xp(-  4a2)  dal  dr. 


t 


dx 


where 


and 


P(e)  =  1  -  P(C) 


Multitone  Jamming 

Since  r(t)  is  the  received  signal,  in  this  jamming  case  it  takes 
on  the  following  components: 

r(t)  =  Sm+l(t)  +  n(t^  + 

where  v(t)  is  the  composite  multitone  jamming  signal. 


N 

v(t)  =  E  /sJ  Sin(a)..t  +  0)/N 
i=l  0 


w. .  is  the 
J1 

Let: 


ith  jamming  tone,  0  is  the  random  phase  which  is  assumed  to  be  0. 


m+1 


(t)  dt 


t^JQ  Sinfuj^t) 


^  Sin((m+I)2irt/T)  dt 


l 

N 


N 

E 

i  =  l 


JTnf 

°  J  o 


SinCw^t) 


Sin( (m+l)2nt/T)  dt 
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The  argument  of  the  summation  is  the  same  expression  as  in  the  CW 
jamming  expression  with  a  single  tone.  Let 


v  .  =  /T7T  /  Sin(aj. .  t)  Sin ( (m+1)  2rrt/T)  dt 
mi  o  7  n  Ji 


1 


N 


"V  ”  N  .  Vmi 
i=l 


The  probability  distribution  function  for  rm  given  is  true  is: 


f  ,  x  =  N(if  p,  ,  +  J  ,  ) 

r  H.  k,m+l  Jm  o 

m1  k 


where  j  is  the  composite  average  of  the  N  jamming  tones.  This  distribution 
has  the  same  form  as  the  single  tone  distribution  in  the  CW  jamming. 

Probability  of  Bit  Error,  P,(e) 

For  comparison  purposes,  it  is  necessary  to  normalize  the  probability 
of  symbol  error,  P(e),  to  its  "per- information-bit"  error,  P^Ce),  (17:711). 


Pb(e)  =  y  P (e) 


This  relation  holds  strictly  for  orthogonal  signals.  The  P^CC)  follows: 


Pb(C)  =  1  -  y  P(e) 


=  P(C)  +  y  P(e) 
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APPENDIX  B:  BIT  ERROR  PROBABILITY  BOUNDS 


This  Appendix  determines  upper  bounds  on  the  bit  error  probability 
for  each  of  the  encoder-decoder  pairs.  For  the  CC1  coding  the  discussion 
is  straight  forward  and  follows  Viterbi's  discussion  (28:239).  However, 
for  the  dual-three  coder,  a  much  more  complicated  derivation  is  necessary 
With  approximations  and  published  data,  an  upper  bound  on  the  bit  error 
probability  is  obtained 

CC1  Encoder-Decoder 

This  encoder  is  characterized  in  Figure  2.  Its  constraint  length, 

K,  is  3,  and  its  rate  is  With  no  loss  of  generality  and  as  an  aid 
in  simplicity,  the  input  source  bits  to  the  encoder  are  all  zeroes.  A 
bit  error  now  represents  a  "1"  at  the  decoder  output.  Figure  B-l 
illustrates  the  trellis  diagram  for  the  CC1  coder  showing  the  Hamming 
distances  for  each  path  when  compared  with  the  all- zeroes  path.  The 
decoding  process  creates  an  error  event  when  it  selects  a  branch  that 
diverges  from  the  correct  path,  i.e.,  the  all-zeroes  path,  at  node  n-m 
and  remerges  with  the  correct  path  at  some  later  node  n,  where  m  is  an 
integer  with  some  minimum  value  greater  than  1  and  some  maximum  value 
of  n.  This  minimum  is  based  upon  the  characteristics  of  the  trellis. 

From  Figure  B-l  at  node  n  a  remergence  of  an  error  event,  from  state  b, 
with  the  all-zeroes  path  has  a  path  three  branches  long  with  Hamming 
distance  5  (i.e.,  state  a  at  node  n-3,  c  at  n-2,  b  at  n-1,  and  a  at 
node  n) ;  also  at  node  n  there  are  two  error  events  four  and  five  branches 
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NODES  n-S  n-4  n- 3  n-2  n-1 


FIGURE  B-l:  CCI  TRELLIS  DIAGRAM 

showing  path  Hamming  distances  in  relation 
to  the  all -zeroes  path. 


a 


D 

•  \ 


FIGURE  B-2:  MODIFIED  CCI  STATE  DIAGRAM 


This  expression  says  that  there  is  one  path  with  distance  5,  two  paths 


with  distance  6,  etc.  This  result  agrees  with  the  initial  observation 
of  error  events  terminating  at  node  n. 


Now  that  the  events  can  be  characterized,  it  is  important  to  identify 
how  many  bit  errors  are  generated  by  each  error  event.  In  a  similar 
manner,  the  term  I  accomplishes  this  in  Figure  B-3.  This  Figure  is  the 
same  as  Figure  B-2  except  along  each  dotted  line  an  I  term  is  added. 

As  the  decoder  simulates  the  actions  of  the  encoder  and  follows  one  of 
these  paths,  it  will  output  a  "1".  This  is  a  bit  error  for  the  all-zeroes 
input  to  the  encoder.  Therefore,  as  the  error  event  traverses  Figure  B-3, 
the  accumulation  of  I  factors  will  indicate  the  number  of  bit  errors. 
Likewise  the  generating  function  takes  on  the  following  form: 

T(D,1)  =  D5I  +  2DSI 2  +  4D7I 3  +  ...  +  2k‘5DkIk"4  +  ... 

=  D5I  (— - — -) 

1-2DI 

With  T(D, I)  both  the  Hamming  distances  and  the  bit  errors  can  be 
counted.  But  what  would  be  the  probability  of  these  errors  happening? 

To  answer  this  question,  the  bit  error  probability,  P^,  is  defined  as 
the  expected  number  of  bit  errors  in  a  given  sequence  of  received  bits 
normalized  by  the  total  number  of  bits  in  the  sequence.  To  upper  bound 
P^,  the  error  probability  per  node,  P^,  must  first  be  determined. 

Because  of  an  error  in  the  received  data,  the  decoding  process  will 
select  an  error  event  over  the  correct  path  when  at  that  node  n  the  error 
event  has  a  lower  Hamming  distance  than  the  correct  path  from  the  received 
data.  By  employing  the  union  bound  over  all  possible  error  events,  the 
initial  bounding  expression  takes  on  the  following  form: 

P  (j)  <  E  P[AM(x!,x.)  >  0] 

^ ^  J 
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1 

where  i  belongs  to  the  set  of  all  possible  error  events  and 
is  the  difference  between  the  cumulative  Hamming  distance  of  the 
correct  path  and  that  of  the  error  event  over  the  initial  unmerged 
segment  of  the  error  event.  But  this  P[AM(xj,Xj)  >  0]  is  simply  a 
pairwise  probability,  P  ,  for  two  code  vectors  over  the  unmerged 
segment.  Letting  the  Hamming  distance,  d(xj,x^.),  equal  d  for  these 
two  code  vectors,  Viterbi  states  that  P^  can  be  bounded  by  the 
Bhattacharyya  bound  in  the  following  inequality  (28:244): 

Pj  <  exp  [  d  I  /p(y|input  =  "l")p(y | input  =  "0")  ] 

<  zd 

where  Z  =  Z  / p(y  |  l)p(y|0) 

Y 

Let  the  Hamming  distance,  d,  of  the  incorrect  paths  be  designated 
as  H(d) .  Then  the  error  probability  per  node  can  be  bounded  by  the  sum 
of  all  error  events  with  all  Hamming  distances  times  their  pairwise  error 
probability,  i.e.: 

OO  00 

P  <  Z  H(d)P.  =  z  H (d)  zd 
"  d=df  d  d=df 

where  d^  is  the  shortest  error  event  Hamming  distance.  This  summation 
looks  very  similar  in  form  to  our  generating  function  T(D,I)  when  1=1. 

00  J  r  d 

T(D)  =  Z  2D 
d=df 

Therefore,  to  use  T(D,I)  as  the  bound  for  Pn,  D  must  equal  Z  and  I  must 
equal  1,  since  11(d)  already  equals  2d  The  final  bounding  expression 
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D=Z 


Pn  <  T(D) 


E 


.a- 


d=d 


f 


This  holds  for  finite  sequences  of  input  data  as  well. 

By  weighting  each  term  of  the  union  bound  by  the  number  of  bit  errors 
that  occur  in  the  error  event,  the  expected  number  of  bit  errors  can  be 
bounded.  This  means  counting  all  the  "l"s  that  emerge  from  the  decoder 
for  each  error  event. 


E  (n  (j)}  <  E  I  i  H(d,i)  P. 

i=l  d=df  d 


E  {n  (j)>  <  E  E  i  H(d,i)  Zd 
i=l  d=df 


But  H(d,i)  is  just  the  number  of  divergent  paths  from  the  correct  path  at 


node  j  and  distance  d  with  i  l's  in  the  decoder  output  over  the  error 


event.  This  corresponds  to  the  derivative  of  T(D,I).  Therefore,  the 


expectation’s  bound  becomes: 


E  (nb(j)} 


1=1, D=Z 


<  I  (d-4)  2d'5Zd 
d=d^ 


where  E  {n^Cj)}  is  bounded  by  the  partial  derivative  of  T(D,I)  with 
respect  to  I.  The  probability  of  bit  error  is  just  the  expected  number 
of  bit  errors  caused  by  an  error  event  initiated  at  j  and  is  bounded  by 
the  same  T(D,I)  function: 


Pb(j)  =  E  (nb(j)}  <  E  (d-4)  2u'3ZU 

d=5 


For  an  AWGN  channel,  as  in  our  model,  Z  equals  exp(-  ^— )  (28:246)  and 

'  o 

the  bound  is  shown  in  Viterbi  and  Omura  to  take  the  following  form  (28:248) 


Pb  < 


/  /  10E  ~  <*>  /  E  \ 

Q  (  /  -fr1) E  (m+1)  2  exp(-mfr) 

/  o  m=0  o 


Dual-Three  Encoder-Decoder 

As  mentioned  before,  this  coder  is  far  more  complicated  than  the  CC1 
coder.  This  will  become  evident  in  the  following  discussion. 

This  encoder  is  depicted  in  Figure  3.  It  has  a  constraint  length  of 
6  and  a  coding  rate  of  3/6  or  h-  Again  with  no  loss  in  generality  and  as 
an  aid  in  simplicity,  the  input  source  bits  to  the  encoder  are  all  zeroes. 
A  bit  error  corresponds  to  a  "1"  at  the  decoder  output.  Figure  B-4 
illustrates  the  trellis  diagram  for  this  dual-three  coder.  The 
accompanying  table  shows  the  corresponding  Hamming  distances  between 
each  transition  path  and  the  all-zeroes  path.  The  table  identifies  each 
transition  by  a  state  at  node  m-1  and  a  state  at  node  m.  Referring  to 
the  trellis  diagram,  at  node  n  a  remergence  of  an  error  event  with  the 
all-zeroes  path  has  7  two-branch  events  with  distances  4  (a-b-a  and 
a-c-a) ,  6  (a-e-a  and  a-f-a),  8  (a-d-a),  and  10  (a-g-a  and  a-h-a) . 

Immediately  the  complexity  is  evident.  There  are  49  three-branch  events 
with  varying  distances.  The  modifying  of  the  state  diagram  in  Figure  3c 
should  lead  to  a  closed-form  expression  for  all  possible  error  events 
terminating  at  node  n.  Figure  B-5  is  this  modified  state  diagram.  The 
states  are  the  vertical  lines,  and  the  lateral  line  segments  are  the 
transition  paths  of  the  state  diagram.  The  terms  D  and  I  have  the  same 
meaning  as  before.  The  product  of  the  D's  and  I's  closest  to  the  end  of 


FIGURE  B-4;  DUAL-THREE  TRELLIS  DIAGRAM 
and  Hamming  distance  table 


FIGURE  B-5 :  MODIFIED  DUAL-TUREF.  STATE  DIAGRAM 


each  line  segment  (or  transition  path)  are  the  transition  parameters 
associated  with  the  transition  from  the  state  on  the  opposite  end  of 
the  line  segment.  As  an  example,  the  transition  from  state  d  to 
state  e  has  the  D3I  product  associated  with  it,  while  the  transition 
in  the  opposite  direction,  e  to  d,  carries  the  D4I2  product.  The  set 
of  simultaneous  equations  that  will  lead  to  the  closed-form  expression 
are  in  Figure  B-6.  Because  of  its  complexity,  the  function  T(D,I) 
was  found  for  only  the  error  events  having  branch  lengths  of  two  or 
three.  This  approximation  is  reasonable  to  understand  the  development 
of  the  bit  error  probability. 

T(D,I)  =  2D4I  +  D5I  +  3D6 12  +  4D7I2  +  D8(3I2  +  41 3  +  I4) 

+  D9 (31 2  +  SI3)  +  D 1 0  (41 3  +  41 4  +  21 5) 

+  D1 1  (71 3  +  SI4  +  I6)  +  D1Z(I3  +  I5  +  I6) 

+  D13(I4  +  I5)  +  D14(I4  +  I5) 

The  general  form  (28:246)  of  P,  for  a  rate  —  where  b  is  the  number  of 

bn 

input  bits  to  the  encoder  and  n  is  the  number  of  bits  in  the  codeword  is 


PbCj) 


FE{nb(j)}  *  F 


3T(D,I) 

31 


1=1, D=2 


E 

For  the  AWGN  channel  Z  =  exp(-  )  .  Therefore,  the  bit  error  probability 

o 

for  those  error  events  of  two  or  three  branches  is  bounded  in  the  following 
inequality: 


+  3  exp 


(- 


14E. 


o 


) 
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The  general  form  including  all  error  events  takes  the  following  form 
(28:248) : 


o 


This  Appendix  has  looked  at  both  encoder-decoder  functions  and 
has  arrived  at  a  bound  for  the  bit  error  probability,  a  measure  of 
performance  in  a  convolutional  coding  scheme. 
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APPENDIX  C:  SIMULATION  PROGRAM 


=  ’jr-  AM  REM  EL 


PROGRAM  V1MA?LES 

DC  »  DUTY  CYCLE 

CNT,UJ  t  G  C.  3  K  <  E  Z  P I N  G  CDJ  *TERS 

r  «  j  ;  Pco<<ErPirJG  cdjuters 

if  c  t  cr jvolutioal  encoded  identifier 

I l!  SEX  l  r£MO2ULPT0r  SIGNAL  IDENTIFIER 

is  ■»  t  number  op  iod  git  b.ddks 
J'  t  IrENTIcIEl  ro  .  jammer  type 
jr  J  JAiUE?:  POWEh 

<  »  NUMBER  of  SIGNALS  SENT  PEP  CAT  £  BLDC1' 

M  t  ENCODER  BIT  MLHOkY  SI7S 

N"  *  NOI^R  °OWER 

N'J  *•  I  MIM^E*  OF  CW  TONES 

S°1  f  EE£0  FO=  THE  AWGN  GE  NE  RA  TOE 

5".?  I  S r  ED  FOR  A  WANDON  Jo  “’HER  NOTES  GEM  "  -  A  70  r>  PQ  )M 

3-D’*  I  SEED  FO f  A  KANDON  JAMmE-v  NOTE:  '‘-t  V  r '  A  TO  r,  SB  J ^ 

S  nu  I  St£D  FO.'.  A  RANDOM  PHASE  GENE  *  PT  Ur. 

T  t  TIME  INTERVAL  r0R  INF  EGRAT I  ON 

WN  i  A  WON  k.’.NOOM  VARIABLE 

DOUBLE  PRECISION  SDi ♦  SD2,  S03,  S3* 

I  NT  EG  ER  “,I, J, I  NO  EX, K,ICC,NJH, JAMMER, CNT 
REAL  NC,DC,J.  ,T 

PROGRAM  ARRAYS 

D  t  OUTPUT  PIT  GTRr.G  3ER  STATE  (‘■'AX  SI^E  51 

C  l  CG0E3  DATA  ARP  A  Y  (it*  Sin  2E>) 

CM  t  RECEIVED  C  0  j  E n  MESSAGE  (MAX  SITE  ?=M 

CC  M  t  CO  N  V  0  L'  i  T  *CNAL  CODE  STATE  T  A  D  L  t  (M^  X  Sin  ?.,) 

IN  DATA  i  INPUT  DATA  (MAX  SITE  120 

JN  t  Ja.MMZk  NOISE  PSNDOH  PROCESS  (-*X  SI7:  :) 

R  t  RECEIVED  SIGNA.  ORTHOGONAL  CO^PORT'rS  (MAX  Sin  ?) 

RH  C  *  CORRELATION  MATRIX  (MAX  SIZE  nXC> 

S  »  T -  A NS MIT TED  MODULATED  SIGNALS  (MAY  SI7E  1;  ) 

wj  t  jamming  frequencies  max  size  20 

OHARA  CTER  9  < ..  1 6)  «  m ,  CC  M  ( •„  1 6  T)  -  £,  CM  (  C  *  2?  t  )  *  1 

INTEGER  Rmq  (■  ti,  U)  ,  1N0  ATA  ( -3  *  12 1 )  ,  ),  S(i  »i"  ) 

REAL  R  ( •_  1  b  )  »  JN  (  *  /  >  »  WJ  (  u  t  ?,r  ) 

»»  **+»»*♦  %*#».»«.^*.»»»*k**4****-»*««  »»*♦*  *  -  «  *  ♦  •  ♦ 

«  t  I  I  »»»»»#»  i  ^  « 

INITIALISATION  OF  PARAMETERS 

THE  FOLLOWING  DATA  ST  AT  PHENES  INITIALIZE  TMf  J  M3DR*  i'JT 
AR  FAYS  TC  ALLOW  PROCESSING 

DATA  (  INC  I)  ,i  =  .  ,7)  /*  "w/ 

DATA  (C  (  7  )  »  3  =  „  j  2  ">  )/?cl 

DATA  (CC!  C  )  ,J*j,  S3)  /to**  *  *' 


116 


r-ATi  (i»»:  ATf  (  1  )  *  1  «-3 , 12.‘  )  »  (  (  -'H?  (  I  ,  J)  ,1s  ,  *  )  ,  J  =  ,  a  )  /  I  2 
•T.TS  (>rO  (’  ,1)  ,1=  :  ,'M  /  9*1/ 

^'Tr:  <*,*:> 

no«(»T  (*1‘) 

IHT7ULI7E  TH-_  RANDOM  NUM3ER  7 T ^ i*  1  0«'  s 


»»****•**;*■***  »  4‘  *■*****-.,  *  .  (  -.  I 


'  L  T  H"  Tnoi.it  VARIABLES  1 


tr  •»  » 
l"  r  t 

ja  : 

v  ? 

J 

r  * 

D“  V 
N  ’  I  ►  r 
WJ  i 


N"i=iK.  OF  i:T  r  I  t  3  t  f  A  5LCC<s 
CC  4  J  Du  )T  IC'NAL  2W)i\ 

Jfr-I.  TYPE  IDE  Ml  F*  f  ^ 
b-T-IZZ  3  CHER 

JAMM£r.  POW~;< 

TIME  INTERVAL 
TUTY  CYCLE 

M'lBsn.  OF  JAMMING  f  3  S  E  S 
Jtl-ING  FRLOUEMCY  A-iAY 


,  cND=?  ’  3)  IST,ICC,  JA  NT ,  JE  ,T,0:,^UM,  CWJtl  >,1*1,  MM) 


*  «  *»"»'»•  1.114,14 


*****■“***•«-*■»*'•  14tM|M»>imtl411>i'<«lllM4lil>l 


PREFACE  th£  -ANOCM  T'i^JT  DATA  STUt.r. 

TM  r  A  T'  fTOPES  TM-:  3L0CK  Ze  0  AT»  ^^ICH  AM  AIL  7£-.OEC, 

r r  a  r:,;  jccu.;s  mT  the  dut3ijt  ah  erfda  is  o? ^ ect ~ r. 


ecuc cMt c y  THE  input  data 

w  *:r-  c  >  is7,i„c,7 ,nd,j*  r,  ji  ,  ci,  >  *j* »  (wju>  ,t=i,m'j  o 

-D7M*T(/7i  ,"I  ►PUT  DATA  BLOCK  HA,  l.  ELEMENtC,  ALl  ^  E r  T  £  >  i 
Tl  ARE  ",I3,"  INPUT  DAT  A  3.DD<E.”A 

ri* ," rv'\io*.\m cncl  coder  nuk-mv,  i»  ,*•  i>  siouestf*," / 

iHE-  1/fL  I S  *' ,  F 1  „  •  E  , "  DiDDNDS.’V 
7  1  T  ,  ••  ft'JC-N  rows-  PARAMETER  s",Fi:,-/ 
ri  .,*•  JAR-EP  hU*'3c3  ”,14, •'  IS  nENl  IFIEO."< 
rt-,"  JAM'S  FOR’t-  FAHAR-T-P  f  ••,-15. 5/ 

r i  • , “  r»n  y  cycle  : M , f  a. . o  / 

F. C'A  DC  ACTS  **  ,  I  * ,  •*  c»  5  1 C  NA  L  S  ,  **  / 

'i  i  i  r- '  \  l  3  a  r  :  "/,T?  .1  :1.  .’//) 
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:o^V0LUTI  Cfv*.  L  ENCODING 
ICC  IDENTIFIES  WHICH  CO^EE  13  T 0  BE  USE*- 

1  «  CONVOLUTIONAL  ENCODER  ,  CCl,  F‘CM  ^CELIE  C" 

2  I  CONVOLUTIONAL  ENCODER  -ROM  T  -IE  fire  DOCUM«*NT 

1 

TMC  c  A  R  r  A  Y  CONTAINS  TH1"  C  0  ">  £  D  31  T  At-.?  AY 

Ic  (I CC.EO. 3 )  THEN 

CALL  cKCQRUi;  ,1  N  D  AT  A  ,  C  ,  Z  Z  1,  M) 

ELS? 

CALL  E>cr^2(l..  ,INCAT»,C,~Cf  ,M> 

ENOIP 


»  #  • 


i  #  »  ■  i  .  •  «  «  a  i 


4M  ’  ‘  -••*♦** 


3-fFY  A'SK  MOO'JLATOR 

3  ?  STORES  EACH  T kAKSMITTE  O  SIGNAL! 

THE  VALUE  IN  S  IDENTIFIES  WHICH  =11  STRING  WAS  SENT 

<  a  ?  ' 

OD  2r  C  J=l, K 

2C  S(J)  =  C(  3*  J-2)  *U  ♦  C(3'J-1>»2  ♦  C(?'J) 


*  *  •  **»***<#*-i.#»»*«***'*-.-*»-+*  **».*.,•*«**.***:«  «•»»**»»*  ‘‘  '  *•’ 


JAMMERS 

JAMMER  IDENTIFIES  Tw  E  TY-E  OF  JAMMING 

1  1  Ei  DACE*  ND  JAMMING 

2  !  SWITCHED  BROADBAND  JAMMING 
T  t  CONTINUOUS  WAVE  JLMMHG 

L  *  MJLTITOhE  JAMMING 
OTHER  t  NO  JAMMING 


IF  (  J  AHtitR.  E0«  A  )  THEN 

CALL  NWJM  (£0*- ,NUM,  j:  ,T,WJ,  JN> 
FLSETF  (JAM*  CT.  EC.  3)  THlN 

CALI  OWJM (SO- ,J  ;,WJ(l),r, JM) 
ELSEI  F  (  JAr*£R.  EO.l)  THEN 

CALL  S  r  J  • :  (SQ3,j;  f  DC »  J  N) 

ELREI  F  (  JA*-.1  EFC.EO.i)  THEN 

CALt  f°JM(S02,J1  ,  JN) 

ENOTF 


n 

m  9  *  4 1  «  #  *  #  *  *  *  *»•«««  *  *  •  *44*  *«4M  *444.'!  «  »  *  «  4  *  9  *  *  •  *  i  «  4  t  xr  4  r 
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1 


COMPUTE  THE  RECEIVED  SIGNAL  OfT  HC.GOf.PL  C  0  -IP  Cl  ENT  1 

V  STC:  EE  THE  VAL'JE  CF  TO  Awy<  D0UPCNEN7S  G1  J  r  F  AT  ED  rr  •>  < 

THE  RANDOM  U'JMr E*  GE  .TRATDR  SGNDr 

THf  SU5ECFI-T  OF  THE  GREATEST  mJE  0s  THE  5  tREA  Y  ID**1TI~IES 
THE  HOST  LIVELY  SIGNAL  SENT 

:  I  =  1 

?“.*  R  <  ? )  =  -99?  • 

INDEX  =a 
On  -  J=f  ,7 

WM  =  <GGOF(SLl)*<  (N./2)*‘  j.:>  ) 

K(J)=EHG(  J,  S  (I)  )  ♦  WIJ  ♦  J‘J(  J) 

IF  (F(lNOfX)  *LT • R (  J)  )  TH-< 

I N  C  E  Y  =  J 
SNDIF 

4*  CONTI  Vl£ 


*  CEhcOUL AT  OR 

*  THE  DEVCTUL S TC>"  TRANSFORMS  THE  RECEIVED  SIGfU  INDEX  I1~D 

*  THE  APPAC  X]  "AT  EC  BIT  STRING  40  FILLS  1  H  £  ru  A  r  "  A  Y 

TF{iK’r£X«GT«3)  THEN 
CH(7*l-2)  **1» 

INDEX  -It  GEY  -  4 

EL5r 

CH(3»I-2)«*.  * 

ENDI F 

IF(p.S"X.G1  .1)  then 

CH(?«*I  -  1)  =•!  • 

INDEX  =  INDEX  -  2 
EL  SB 

=  •)• 

EOT  c 

IF  (INDEX. GT.  )  THEN 
CH(T«I)r«i» 

EL  Sr 

CH  ( 3*  I )  =fu» 

END  I  e 

* 

*  HE  FE  THE  DI  MODULO  Or  RET'JRN;  FDR  lHr  NEXT  F  EC  El VE  D 

*  SIGNAL  IN  TH£  DATA  '’LOCK 

» 

I  =  t  *  1 
IF  (I  .IE.K)  THEN 
GO  TO  3.1 
F.4  0IF 


»  %  .T  +  .*«  i.  ■.<  i  ,».«  n  «.»»*>  «»*  H  I  »i  >4  .  *..  «  n*.  (  i  i  <. 

* 
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VITFF°1  0 £CQ HER 


rwe  C£LL  OF  T  •-'  £  DECODER  MASSES  THE  NEIES  SAPY  I  NrO  rMATI  0  >4 
TO  FOLLOW  i  HZ  SHJMEST  HAr-MiNS  DISTANCE  THROUGH  THE  :mODcR 
ST ATE  TF ELLIS 

-  c  (IfC»En.l)  THEN 

o' ll  s/n::i ten, ?c- ,cc-:,a)  • 

ELSE 

CALL  V]T0C2(CHf  £lf  ,CCT,^) 

-no  if 

OD  **  J  =  t  ,1!  5 

]F  (0( :  )  ( J«  J)  >EO.»i  •)  F  H EN 
CM  =  C  NT  +1 

END  I  F 

51  CONTINUE 

#  *  f  #•  ***-«'  ^  +  +k+»4  +  +4  I  if.  H  4  4*  »  4S  «f  «  f  #  4A  A«#  «■«  t 


*  3R  1ST  THE  RECEIVED  "IT  SI  TINS  ?ro?EO  IK  *  C  ) 

» 

w^ite  (* » y .  a  . )  CNTji  <  j)  ( s * i : . ) 

O'i  c051*T(/fI-,TC»A) 


Pr  TIIPK  FOR  ANOTHER  DATA  CL0~< 
CONTINUE  WITH  Tm_  PRoGRA- 


IS?  = 

IS7  - 

1 

I F  (T 

S-'.CT  . 

i  ) 

ELSE 

TO  TC 

1. 

ENOTc 

CO  TO 

r  • 

*  #  »  ***■***••«  **»-4**»*«***-*»#*.4-***>***.».*|  f»,4)  •  |  4  I  4  *  *  *  ,1  '  <  -  »  «*..  4  -  •.  «  •  »  *  -  4 

«  *  ■»  *#**4*4  »-»»»***4*««-»-**'*-**4  *•#■****•**>  *  »  •  •  »  *  4  *  *  *  *  *  *  *  #•*■**  *»•.***»»***•*** 

*  EMC  OF  INPUT  OAT; -  SIGNING  OFF  • 

6?  W1ITE  (•  »  D  ■  r  i ) 

915  FO?**!  T(//,Tlt  ,"ENO  OF  INPUT  FILES"* 

STOP 

*  * 

*  *  *  *******  «*»i  4  *  •  •*•#**##*«•*•*•*  -*•»*»  »*«*»-»«  i  >  4  <»>  *>4*1.  «»♦»  l  .  i  *  »  a  *  *  * 

ENO 
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c'»-0  LNCORl  (I, IA,C»3?y» "> 

.  «-»j  a  •  ■>  .  V  4  + 4>**4*«*  *<-«.» 


*■*♦*♦*-*»*»»»****»»»->♦■*♦**»»** 


HIS  :U?KOUTI^£  IS  4  SIHMr  1 0  M  /  0.  UT 1  D  NAl  CC"-'?  w'TM 
a  r.0i~i'~r  Z IKT  LENGTH  =  3  L  40  a  IDDInG  RATE  OF  4  HAH 

TMf  CALL  STATEMENT  PARAPET £RS  ARE 
2  t  INPUT  8LG;<  length 

ia  t  input  3ir  array  (max  SI7E  i z.)\ 

:  i  C‘ >T PUT  3rr  A: vRSY  (MAX  31  ?r  2'  *. ) 

EC*'  t  C JuE  MAT -21  X  HIM  OJTP'JT  C OCpWjir  p  — 

STATE  TRANSITION 
M  J  NUMBER  Oc  MEMORY  PITS 

TNTEGRT  I»M 

tntegef  it(-’:i2r)  ,cr?  23  ) 

CHA-  AT  TEC  CCW  (  ?  A  ? )  *  8 

j-uc  j  f  Ti r  N4  L  PARAMETER,  J,  13  A  SOUnTInG  VA-.I^LE  * 

TUT"'  -  2  J 

H:  INITIALISATION  OF  THE  Coon;  matrix  * 

L  V-  P  Y  £l£mc.’«T  RE  PR'S  RUT  3  C4E  STATE  tnAK3It10N  * 

PI-SI  rl!  IS  THE  PRESENT  STATE 

sEICKD  SIT  IS  the  u-VT  state 

L  A  S7  SIX  3ITS  REPRESENT  HE  ENCJPi-..  OUTPUT  PER  *■ 

T A AO SIT  3  ON  * 

TSMC  )  =  *.  i„l'  1  3* 

"'“'*(  I  )  -  •  2l:  3  tf  '  3  * 

cmi? )  =  *i( k:  i : •  i* 

S'-M"  1  =  *12.  I  J  • 

Z MU  T  =  MU  T .  1  * 

:*?'*(.)  =  *22;.-  i  »  * 

”■<(■  >  =  '21*.  l,.  ’» 

•TU’  >  =  *33*.  „*  1* 

T'^.WS-  u  Gft.E-ATTON  IS  ?  “•OOJ.3-2  OPERA  T I  ON.  PO"  r  A  P  H 

INPUT  CI "  AW  THE  TO  PRE/rDJS  INPUTS  A  HO-  PIT  * 

C  Or  3  WC (  J  IS  FORMED  1*4  THE  rDaOHU&  MA*'4T‘r,  * 

00  1  ■ 

:  <  c  J-l)  -  *'CP  ((  Tf  {  J+i)  +1 A  (J-l)  ),  ?) 

r<2»J)  =  MODCAIJ  +  l)  +  1  A  (  J)  *  I A ( J-i )  1 2) 

CONT I KUE 
“  =  ? 

• 

SUP’OUTrt-E  C  DERATION  OOUEf  EETJRN  TC  MAIN  PROS?  AM  * 

f 

1  *  *  *  *►*■-»<»  •■♦4  4  ♦  *  -  »  »  r  *44*  •  I  ♦  ^  4  4>  »  »  44444  *4  ‘  *  4 


♦ 


CTIr  £  ENCr~2  (i,  IA,  C,COH,  *> 

*  ■*«*''  *  ».  r  *  4  l,  «.  «  «  »  4- A  -*  *■*  *  *  ►  *  ♦  #  *  *  ♦  »  *  •  *  *  •  ■*«» 


T“  IS  EH3*0l'TiN£  IS  A  ;UJ-  E  CO-OLEX  CONv/jl  U7I0  HI 
>JTTW  Z  Oi.E  Tr- e  I  NT  LENGTH  =  AND  4  3  GOT  KG  ivATE  rc 

HEN"  DARAHET "P.5  £*E 
1‘.’ PUT  9LDGK  LENGTH 
1,P'JT  BIT  A  -  -  A  Y  (  H  4  X  S  J  7  E  12  )' 
C'Jl  PUT  3ir  4  .  v  4  Y  MAY  1 1  2:  "r) 

CODE  MATRIX  Wi. T  H  OUTPUT  rCCEWT~r 
STATE  V '?!  US  JT  I  ON 
AMH5"R  0-  <r"IGFY  31 TS 

INT  GIF  1 , 1 

T  JT  :  G-  F  1/  (■  3  !i  2.)  ,Ci  (.  t  2t:  ) 

Tma-ACTLC  CC-(  to?)*8 

T M  S  I'T.F.JL  r  A  r  1  ‘  £T  EPS  A  -'  r  J  <<,  bf'JNTI'G  VFA  =  J. 
iN:E'r"»  j,< 

T u  t  V  IT  lull Z5T1GN  OP  TME  C?OIiG  1 f  TR 1  * 

■  Ji  f  Y  ELE-LvT  P£P3"?EN’i  CNE  STAT[-  T  P  f  t  3  1 T 
r.i.ET  "I"  I T“E  PRESENT  SfAlE 

rE:CN:  3 1 T  is  THE  NEXT  state 

last  Six  JITS  PEP?rEENT  THE  ENCODE-  bUT‘5L)T 
I'  ansition 


•  > 

H( 

) 

= 

*  :  t  -  i  J  •  ;  * 

-  ^ 

*(1 

> 

~ 

•  lu.lrl  • 

•ME 

1 

= 

•. £.i!  1 .  ■  • 

^  *■* 

M  " 

1 

s 

•  JEHU  * 

z-  -* 

-  ( ' 

) 

= 

*  M  1:  .  j  i  1  * 

r  -* 

•( 

1 

r 

•  1 1 . 1 .  1  * 

M 

) 

s 

'•-t  n>  111  * 

"(' 

) 

r 

•  ;  1 1 1 1  •.  1  • 

2  2 

*  (  ‘ 

) 

= 

•it.  1  r  i  * 

~  * 

M  ' 

) 

r 

•ilM  ill' 

t.  -J 

'  (  1 

•) 

r 

•12. 111  .-I* 

r'  ^ 

*  ( 1 

1) 

= 

•i  3.  uni* 

r  0 

•Ml 

2) 

r 

•J* 11  lr  1:  • 

r  0 

Ml 

?) 

- 

•  5 4 1 '  a ■  >:• 

•Mi 

O 

= 

Mtiimr* 

r  ~ 

Ml 

r) 

r 

•1/1131 jf  * 

*M1 

n 

= 

•  ?  1  L  1 U  '  1^  * 

Mi 

*•) 

= 

•21'  1 1\  * 

■»  ^ 

V  J 

Ml 

A) 

r 

•22i  .  -Ill  * 

Ml 

p> 

r 

*231.11  31* 

r  ~ 

M  ’ 

*> 

r 

*  2  *:  1 1 J  ,  31* 

r  ^ 

M? 

i  > 

r 

*2.111'  11* 

**  ** 

“  (  ? 

?) 

= 

*  Z  0 1 1  J 1 3 1  * 

^  (  r> 

2) 

r 

*271,  1111* 

■*  ^ 

M? 

n 

z 

*3n'  11.  11* 

r  •  -“H 

'  ( 7 

-) 

- 

*:•  it  1.' '  1* 

*"*  ** 

■  f  -> 

f ) 

- 

»::•■.  1 1 1 1  • 

- 

•M  ? 

*  > 

r 

•33L  :  Jl'.l* 

•HE  ( 


ALL  STATE 


A  : 

c 

I  Ch  I 


»  *»•-** 

rr  ■(> 

4  h;i_c 


o  r- 


"  l:  g 


o  b 


c  cr 


122 


9 


CCM(2  f> 
ocm<?<3) 
C  C  K  3  n 

COM(  3  1) 


COM(t  2) 

s 

•f  i  1  JlJL' 

C0l<(  7  2) 

- 

•*  ill  nil  • 

C  C  *  (  3  h) 

s 

'A  211  J..  -I* 

CC-*(2  p) 

s 

•*■3111'  K* 

C  “  -*  (  7  f  ) 

= 

Jill' 

"‘(3  7) 

= 

•*  S'  '  11  '  1* 

“0J<  3  6) 

= 

•t;(ii  n* 

:o 4(3  s) 

r 

*  4  7 1 1 1  i  ?!• 

:-“(  r > 

= 

•L-.  If  11  >.  1* 

coun*  i) 

s 

•:  111  51 11* 

cvir-  ?) 

r 

•l  2111  31* 

C0*(4  2) 

= 

•r  3iij  on* 

CCn-,  4) 

s 

* f.  tluilll* 

CC-‘(4  *) 

= 

•  •  5 .  r  j  l  J  l  • 

C  (  .  <■) 

= 

•r&cn  ii * 

C0M(4  7) 

r 

*T  7tl  J,  V* 

P) 

- 

*c  jinn. ' 

r:«(4q) 

s 

•t  mil  ji  • 

OOMC7  f ) 

r 

»f  21.  ...1C  ' 

00^(2  1) 

= 

•f  31 r i:  ;c » 

CC-C  2) 

= 

•U  .  m  i® 

CC  M  ( E 2) 

z 

*rit  i  ini* 

c:m(p  4) 

z 

:i* 

CC*!"  r) 

z 

•ETcriou* 

CCM(r  FJ 

r 

*71111111* 

CCM(?7) 

= 

*7  111  3 1 r  1® 

CCM(=  ?) 

r 

•V  21. 1L  11* 

C  3  M  C  C) 

= 

*731'.  I .  'I* 

C0*(*  f) 

z 

*74(111.:  * 

CCM(G  1) 

z 

•71! lillt  * 

conn  p) 

z 

*7fm: 

CC*MS  3) 

z 

*77 L  on-.  * 

?A111.  J, 

2*  li.u  ll • 
?ti-  11H* 
271.  Jlli  • 


sor-wc-k?  GFNEkATIOIm  is  a  SET  C?  MjOULD-2  equations.  St<  "TTS 
(3  PRESENT  INPUT  A*0  7  rim  GEN  If  ATI  A  SIX-3IT 

OODLWUt-J  IN  THE  FOLLOWING  ASSIGNMENT  ST  A  TS  HE  f-TS  . 


oo  2:  n  j  =  1, 

K  a  K  ♦  1 
C(fc*  K-l ) 
CAb"  <-i  ) 
C(fc‘ K-3) 
CCCf  K-2) 
C(G*K-1) 

CONTI  M3E 
M  =  3 


1*3,3 

MODI IA< J  +  2)*Ifl( J- 1 1  ,2) 

MCG  ( I A  (  J  ♦  i  >*I  A  C  i-c)  ,  2) 
f * C L' ( I  A  ( J  >♦  Jfl  ( J-  *  >  ,  Z> 
KOD(IA( J*l)+iA< J-J  )  ,2) 
HQIKIAT  J*2>*IA(J)  *IA(  J-2)  ,2) 
M00< iA( J  +  2)*IA{ J- 3 >  ,2) 


SU^.OUT!  I  E  0°E  =  ATiQ':  OO'-'L!  \ZTJIN  TC  *'AT'J  do AH 


J  .U  «  ^  4  i  4  I  4  -r  - 


*  .**.*.-  ^  -  H'  U  *  M  1  ^  ^  W  *  *  *  *  4.  ~  -  A  » 


S'l^piiTU  E  (  302,  J.,  JN) 

*  a  »  *  *  «**(*.«»♦***».*.*+*»■»#»**••»■*****»***•»<**  ***'*»»*J1J,*‘  ♦”■**-** 

THIS  SUEfCUTINE  IS  THE  GrOS^C  SD  J  ^  L  F  WITH  THE  FOl.LO-JIJ 
Of  LL  STaTE^NT  P-.r^H£T£3S  * 

ED?  t  NjEHAL  RANDOH  *JJ> 13E\  GENERATOR  S"r 

J.  1  JAHHER  p CM  EH  °A^HETE9 

JN  S  J->  HKI  MG  EDI  Frill  ENT  HR- AY  * 

COUPLE  P'ECISIDK  S1'1? 
t  EAl  J  ,  JM  .  »>  ) 

Twc  iMEFNfl  COUNTING  VA?Ii=*[.r  jS  j. 

INTEGER  I 

THE  JAmnING  CO  EF-ICI  ENl  ArCAY  IS  ="1LLEP  WITH  EIG^*  ML  JET 
WITH  THE  ►'’.AN  OF  AND  VirlONO;  Dr  (J  /E)v*  .  E 

00  5  "  I  = 

„N  ( I )  =  GGNOF  (CD?)  *  (J  /  2)  M  >.r 
CONTT  F’JE 

tME  S'-lofrUTITE  0-  ERA  T  ION'  IS  DOM_f  ki.T'J\N  TO  'HE  MAIN  PRO  I 


*  **  n  *  »  *  ««  *.*«.* 


»  |  I  »  *■##  «<«*  4«  .  * 


B  J  (  S  C  2 , 


M,  a *  »  4  *♦  »»  *»» 


;u-?uu’IN:  IS  i  HE  S^I’SHED  «u^?tr,[  J*-  <*r'-  W*‘M  T4' 
OLLO^Ii'G  CALL  STATcM-NT  PAp.A'IEIERS 

tri  {  KJRRAL  RLtOCH  ^  -  N r  P  A  T  i  s.  t‘--n 

j  j  J„VM£R  Pj*IER  3ARAMc.TcR 

•  DUTY  CYCLE 

j.,  .  r  a  -  ? ;  y 

:  c  -' c L E  Pr  eci sic  n  so* 

r  r ‘ L  J.  , DC, JN C  *7) 

'  -J  :  1 t  T  . *■  K  t  L  °  A RA  *'  E7  £  RS  A  \  £  * 

i  (_  p  *■ '  .*'  ?  T  <  mNS  R I SD I  jM  3i  j  J  <  L  F-- CIRNC* 

T  «  CjUM  ING  VL  RI  a?LE 

~  l  :  CO  US 7  AWT 

;  nANDOM  V*jIA~Li.  /A<.IAS'.,E 


'L  1‘  u  A  ,  P  i,va; 


i :  f  i  l  r  ~ 


ru  f  J fy,;Nr,  COEFFICIENT  *“RAY  J5  FILLED  IN  THE  "OLLO^HG  LOG 


'  I  =  »  » 7 


f!_=  U5  Z  PI  '  li  *  (1*1) 

:h?  variance  i?  calculates  m e xi 


=  J 


(r"  -  ?Iu(i_=Ht  *  r’D)  /  fti  9W:  )  /  ? 


1M-  Jf  mI'-’G  COE  FFIC I  ER"  3  t<£  Oi*LC'J  LMEC  I*  "  ur 
5  7  ~  £ k  E.  N  o  USl'.G  TH"  A’ODY  Ol  NEr  A  T  OR 


■  1C)  =  Gr-:4r  ~  (S  D  3)  *  '■/!•'  *'  '..5 


OPERATION  IS  DONE:  RETURN  TD  THE  HCCM  F’jGRA^ 


»»*.**'>  ♦ 


l  Tit  £  C  WJf  (SC*, Jr, 0MEC5  J,T, JM> 

ru-f  cL»Tifj=  is  the  gdntiyjojs  wav:  jammer  with  the 

r'-'lLrwiNC  CALL  STATEMENT  PARAMETERS  * 

-  O'  t  RANDOM  N’ IMS  EG.  GENERATOR  SEED 

J  t  JE  MMf  *■  SIGNAL  ENERGY 

c««£r,<!j  »  fg.qjensy 

1  }  T.v  A  MS  H I  S~  I  ON  STHE 

JU  f  JAM  M*  J  G  03  -T  Zlz.  NT  A :  <  A  Y 

1 ")  l"LL  Hf  EXIST  ON  SOt 
■? r  •  u  J  ,  0M:GAJ,T,JN<~:7) 

rLj  «  1ST:  ^  NCI  P AKA  METERS  ARE  t 

coy  *  constant  -jn:ti3<  of  p? e cuc.n: i ec 

:  •  CO  UNT  I NG  VARIAR.E 

ChECCK  t  T  ’AM'S  ISSIQN  F?E  3  JENC  Y 
i.  M r  G A 5  t  STORAGE  LOCATION 
~  I  !  CONS'.  ANT 


vE-OL  COM,  °I»  PMEG A“ ,  OMEGAS 
"!  =  ”•  1 '  1:  E  if:  3 

Th  e  Jf.  M  1  s&  COEFFICIENTS  ACE  IMTIALI7:''  T0 

'.0  J"  r  I  = 1,7 
jm  =  r 

rui  c  f'  L  L  0  W  T  M  G  LOJP  FILL?  "HE  JAMMING  C  0  E  F  FTC  T  rMT  TRp.AY 
m  ’i.  *  1  =  1,  3 

0-ECAJ  IS  COMPARED  TO  EACH  TRANSMISSION  F9  Zn  LE  JC  Y 
T^  (QHLGAJ.  3)  THEN 

J„C)  IS  CCLCUl  A  TEO  IF  OHESAJ  EQUALS  Ok'  "Ft  H 
Jt(I-l)  =  ((J'»*T)**:-.5)/2 

CLC  E 

t  jN(i)  is  calcult:*  th:y  a ■>.:  not  ro'«\ 

Of  EGA 3  =  OficC-AJ  *  2  *  =1 
of  eg  a  h  =  i  *  o :  .  *  *  •  =  1 
CLN  =  {2*  JMEGAM) /( OMEGA  5s  *2  -  0*"  EGAH-  ”2) 

Jf(l-l)  =  ( J*  /!)*♦  .••••:D'i*3IN(:H-GAS'T)/E 

'  ACT  F 
C0I|TT  wj" 

G ' »  ~  '  l  H  T I  f  E  OPERATION  IS  D  JN£!  P'TJP"  TO  THE  VI'J  cR0G'<-\ 


juti»  e  fwj*  (*•:..  ,*.u“,  j , r  ,wj,  js> 


TH  7S  SUSFPUTINE  IS  THE  H'JLTITONi  JA^urR  WIT^  T 
rC>  LLOWlf.'C-  CALL  ST  ATE  M  EN’  * 

ST-  :  RAN  DO*  N'JMOE.  GEnExAYQ*  SELC 

N’UM  t  N'J  M3C  R  0^  JA" II  VS  FVFOUE  NS  Iff  * 

J.  J  JAMMING  SIGNAL  S'! S'. GY  ; 

T  !  TRANSMISSION  TIMS 

WJ  t  Jh  VMI  NG  FREOUENOf  A  KF  A  i 

JN  t  JAMMING  COrFFIGISSir  A  c  A  Y 

COUPLE  PRECISION  SSL 
REAL  J  ,T 

\E*l  J*  C  !7)  ,  WJ  ( ;. :  2  .  ) 

INTEGER  N  ’ J  R 


THE  INTERNAL  PAAA^ETLRS  AOS 

I,  Y  t  COUNTING  (/ARIAS.SS 

V  J  INTE<- MEOTATE  jamming  COEFFICIENT  A  r~  j  y 


I NT  EG SR  7  ,* 

REAL  Wt:7> 

THt  JAMMING  COEFFICIENT  AK-.AY  IS  CALCJ.ATEO  tm  fwr  FOLL  ON’  N 
LOC3  Fy  CALLING  THE  CWJ^I  SU-RJJIlNd  A  N  j  ft  VERGIN':  THE 
COEcF:CIFKTS  FOR  EAC  1  JAM  II  NG  frEOUlNOC 

DO  5  0  CO  I  =  1,N<J“ 

CALL  CWJH(SD4, Jl ,WJ(I),T,^> 

CO  ElwV  K  =  7,7 

Jf  (K)  =■  J  J  (K)  ♦  «(<)/ 

CONTINUE 
CONTI  NUF. 

StlE^CUTINi  CASSATION  IS  DUNE?  RETURN  TO  THE  MAIN  f  JG’’.  A  M 


.»■**»*♦*« 


<*♦**♦•*»*►♦♦*%♦->«**¥■•***«*»• 


1  • 


S'l^-UiNL  V ITC3l  (CHjR.CCM,3) 


\ 


4 


t 


-  *  »  »  »  *-.»*»**».*3  ♦»*.**.*♦♦*»,  <***►►*,*4*,  *.***>  >  , 

thi~  su=rci,nui  is  the  vite^i  jeecp-e-  f?*  n-- 

CC  1  rnfJVC  tUTi^NAL  CODER 


'  LL  SUT 

“ M -  NT  PA 

RAH; 

7  "NS 

&  RE 

t 

EH  t 

E3TI*-m 

Tin 

"03£3 

SYH30L 

A  O  AY 

K  t 

1,'IH-t  - 

Qf 

£L- HE 

HTS  IN 

3  H 

ECU  « 

to  ns  itioni 

HA  r  R 

I< 

3  t 

C.FTIHA 

TEn 

sojrc 

;  iN'Ff 

O'ICN  AA  “AY 

riTEGCF  K 

nr 5 AST  £?•  3  ( .  **»  )*i  ,C3'M:  se,’)  *  d,CM( .  \z\ : )  *< 


TM  6  INTERN* L 
OP  t 
HMG  I 
3  t 
I  FH  I 
I7G  t 
CNT,  J, 
I A  i  I F1 


PARAMETERS  ARE 

SCr.ATC.H-'OO  3HARA37rr  e  •  f  tv 

the  h*nmihg  3isiv*:s  a  5 «  A  Y 

(OMUL  A  T I V  E  0  IS  f  ATE;  tR:tY 
P-iFSf  NT  TRANSIT  I  3T  S  T  *  T  £ 

N-"  XT  TnANSiTIOV  5  T a 7  r 
JCH  J  INTEGER  20JTTINC  X6AlAnLT 
TEMPORARY  iTO?£Gi  OC  A  T  I  ON  S 


Chara  CTER  FT  (i.  o  )  *1 :  l 
I  lJT  -G  H  MD  ( l  s  ‘  ,  „  to  ,DC  t  h,  J  *  1  •  3 ) 

INT£CC?  j ,:a ,lr » iCM, 1T0,CNT,ICH 

* 

•  PAK4METEr-  I  NI T  IAL I  ?*  T I  ON’ 

* 

3A'rA  n(-  , u) ,  D(i»  )» 
no  2 ; "n  J  =  .  ,L 

MJ)  =  •  * 

Fo  <  J  >  =  *  « 

23"  COST  I  N'lt 

C‘£T  =  ? 

23  *'  1C«  *  1 ,  ♦ 

?’  e-  J  = 

TO  2  :  £  L  =  ■ , 

HKD(J,L)  = 

3'JILUING  TH£  HAf  **F'G  DISTANCE  ARRAY 
21".  J  =  - »  ? 

]Fn  =  ICHARCOC  « <  J)  (i  ti) )  -15 
JTO  =  I  CHf  R  (CCH  (  J)  (2t2)  ) -15 
if  (ch(igh)  ,  cj>  cj  tr)>  them 

HN  D ( IFm , IT  0)  =  1 

F‘JOIF 

IF  CCH(ICh*1)  »•»  c  •1C-M  ( J)  (FtS)l  THEN 

HM0(IF^,ITO)  =  H-D  (IF-J,  I  TO)  ♦  1 

EWDIF 
l"  CONTINUE 

THE  /  lGOF  IT  ^  PkuC£:  C’G  "v 

-  PiTfR-'INING  THE  MlJi'l'  HAN.-iSG  TO  ST  A  NT 

FOSTinLE  STATE  T  O  NS  j. «  i  3N  F- D"  S'  FT  H  ENT  '3  3 


30 

30 

20r 

30 


3(2,0,  3(3,0  /  i  ,9c-c,?9',pO/ 


2»  2 


7h 

an 


128 


-  SAVING  the;  MNIMU"  OISTANSE  IM  THE  D  P.'.tv 

-  CONCATENATING  T-M  C'J^EiPOUCIVS  1  Cr  '  "M  HE 
F  r<-  Y  AT  EACH  *> T A ” E 


H  =  C(C,CnT)  ♦ 

HMD(k|i ) 

10  =  C(1»CNT)  ♦ 

H  H  C  ( 1 1 1  ) 

IF  (Tfi.LE.IF)  TH 

FN 

r(c  ,CnT*1> 

=  IA 

eU>  =  \* 

//  3P(c) 

~  W  .>  Z 

[  (n  CM  +1) 

=  13 

F  (  C  )  =  • '  • 

//  3 P ( 1 ) 

■‘nif 

I A  =  r{2,CN’)  ♦ 

H-C<2,1) 

13  =  r(3,CNT)  + 

h-D( 2,1 ) 

IF  (I  A.LE.I  c)  TH 

IEN 

CdtC'vT  +1) 

=  IA 

F<1>  =  •  • 

//  3 P  (2  ) 

ELS  c 

C (1 »  CNT ♦! ) 

=  13 

!(1»  =  »  » 

//  0P(?) 

ENOIP 

i a  =  r(f*,cn>  •*■ 

|J-D<  .  ,2) 

13  =  C(l,CM)  ♦ 

HHr-(i,2) 

IF  CiA.LE.IF)  TH 

:  i 

C ( 2  »  CnT  +1 ) 

-  IA 

F ( 2)  = 

//  3P(l  ) 

ELS  E 

"" r<2,CKT*l) 

=•  13 

E(2)  =  *1* 

//  3C  (1 ) 

EMOIp 

IA  =  rtz.Gtn  ♦ 

HM0(2,3) 

13  =  T(3,C\T)  ♦ 

Hjn< 3,5) 

IF  (ia.LE.IF)  TH 

r  N 

n(3,C>;T*l) 

=  I ' 

e  ( 3)  =  *1* 

//  3 P  (2  ) 

ELS  F 

0(3,CNTM) 

=  13 

r  ( 3)  =  *i» 

'/  3P ( 3 ) 

THE  FOLLOWING  LOOP  RESETS  THE  55  ARMY  FOF  THF  “OVE  * 

TO  THE  NEXT  JEOH  NHH9ER  CHT  * 

* 

DO  2r  £-•  J  =  c,t 
3  =P(J)  =  E(J) 

* 

INCREASE  THE  OEcTH  NUH3E*  CNT  * 

* 

CNT  =  CKl  ♦  1 

90  CONTI  HR 

K 

SUF^C'JTiFE  C-P  ERAT 1  ON  IS  TONE:  «STU.\N  [D  TM[  HAIN  F-.OSTH 

* 

*  |  |  »  »  »  *  J  *>#  »  i*  *4«  .  IM  •  *  :  »  .  «  • 

1>'D 


•»  i  ' 


ITI*  £  VI TDC?  (CM  ,  K,  CCM, 3  ) 


r^is  susroviiKZ  is  the  deco de-  pur  r wr 

DUAL  THREE  CO*jVCLUTIOMAl  cpoer 

rMf  GALL  s T /  T  £~»  =*V7  P  A  R AUE~  ER3  ME 

:w  »  E3  TI" A  T  ED  CODE.:  SYMSGL  A  PF  A  Y 

K  s  N'JMRj.^  Or  Icp1isrs  If.  :h 

CCM  J  T  A  NT  IT  lnS  MATRIX 

b  \  lSTIHATEO  SOJRDF  iNPr  Ktt’ICK  4?-  AY 

TNT  EGEF  * 

CM A=  ACT  £f  G  H  ( !  2‘  v  )  *  i » 3  ( •„  »  B  )  »  1 i  5  ,  C C  M  (  .  S  f.  3)  *  5 

THE  INTEFNAL  PARAMETERS  ARE 

DM  »  COO EWOrO  STORAGE  LOCATION 

?P  t  SCRATCH-PAS  C  -<  A  RA  D  T  “  '■»  f  -AY 

JM3  f  T^E  HAMMING  DI3IAS3E  A '  -  A  Y 

D  *  CUMULATIVE  Dism:-  AY-AY 

CN1,  J,  ICH  t  INTEGER  C'JjJTli.G  If  A  fRi  uEL  E  5 
IFF  I  PRESENT  TRANSITION  STATE 

I  Y  L-  t  G'YT  TRADITION  STATE 

X  t  TEMPORARY  STORAGE  OF  C  J”  UL  mTI  7-  DISTA4CE3 

csv  t  present  state  storage  ^oca-yion 

CHARACTER  OF*!,  cp<  tS)*10? 

INTEGER  HMD  (».**,.  J  ft)  ,  D  <  y  t! ,  C  *1, -•> 

INTEGER  J,CNT»:FM,iTO,T,X,:S  V,IH 

PARAHETEF  1KIT  IALl?ATIOk’ 

data  nu,.  >, j=;,a>  /  .,5 *"<-,/ 

OD  ?  ?  J  =  I  ,e 
F(J)  =  •  » 
ro(J,  -  .  # 

CONTI FUE 

cvr  =  ■*» 

O'*  3r  T  ICH  =  1 ,  <,  fe 

CO  33  .  J  =  j,c 

OC  33 j.  1  =  '  ,  8 

HM3<  J,I)  =  3 

CM  s  *  • 

DO ►SITING  THE  CODEWORD  IN  DM 

CO  AD.  I  I  *  3,3 

OM  =  OM(2f  )//CH(ICM  ♦  I> 

CONTINUE 

OMILDING  T Hr  MAF’MiNG  DISTANCE  ARRAY 

co  ^  I*  t  r  —  .  1  u 3 

JFV  :  TCH'  -  <ro-*m  (  i !  1 )  »  -15 
1*0  =  7rH-,  (GCMD  <2?2>  )  -15 
DC  1  ?y .  J  =  1 ,  S 

I F  ( D  •'  (  Jt  J)  •  >r  #  CC 1  (  l  >  (  J  *?  *  J  *  3)  )  H'v, 


>*  » 


HMD  (  I  FM,  I  TO) 


' 'OIF 
CCKTJ-JUE 
COKT INUF 


Hlu  (Irr  ,ITC)  ♦  i 


TH  *  flGOFIT^I  Bf<r'C ETDS  nY 

-  OtTf  ■<"  INING  TWt  •'J.NT-J-  DISTANCES  FOt?  -A? 

F  C  S  f  I q  L  L  STATE  T*AK'3iUj(  Ft  o-  hi  FTH  Crj  T0  CNT+1 
*  S  A  V  i  N  C-  TH  E  ‘,*UIk'UM  DliM'jCL  I  K  'HE  f  A  r  12  Y 


FC‘  CAT  EI.*7I»vG  T^i  CO  H  i  E  3  =*0  NDi  K  C-  1  0:- 
<?  r-  < :  y  at  r  a  c  h  state 


«I  ■*-< 


*7 


CO  «,?.  i  I  = 

csv  -  * 

x  =  r- (  , CM)  ♦  h^t(  ^,d 
DC  v:  J  =  1,7 

I  F  (X.GT  ,G(Jt0t.T>  ♦H-'O  (  J,I)  )  T  WEN 
X  =  OCJ,  OUT)  *  HID  <  J,  :  > 

CSV  =  J 

ENOI*7 

cn.::  j'j: 

j  =  CSV 

P  <  I , r 7  +1 )  =  X 


IF  <  r  S  V  .  r,j  .  3 )  T  H 

irv 

0(1)  =  *1* 

// 

tJO(  J) 

CSV  =  CSV  - 

4 

ELSE 

9(1)  =  •]* 

n 

9=(  J) 

FKCI' 

IF  (CSV.GT.l)  TH 

EV 

9(1)  =  *1* 

// 

B(I) 

CSV  =  CSV  - 

ELSE 

3(1)  = 

n 

’(i) 

=  ».T  IF 

IF  (  C  S  V  •  GT  .  ’ )  TH 

~~  ii 
% 

3(1)  =  'i* 

// 

5(1) 

ELSt 

9(1)  =  •' * 

// 

9(1) 

Ef  DIF 

COM  If.l'E 


THE  FL-LLCWlrO  COjP  ~  ES£T  5  THE  ?=>  AkkAY  FQk  THF 
TO  THE  NlXI  DEPTH  NUhOc’  CNT 


h  rv“ 


no  •"  J  =  c  ,  7 
S3  ( J )  =  P(J> 


INCREASE  IMF  D  c  P  T  H  MJM9C*  CM 

CM  =  CNT  ♦  1 
COMT  KHF 


S't*-?Ql»TIh£  r3rcATIOI.  IS  ?jv:t  VET'JPI.  TO  THE  HtlM  ( 

f  ■«>.,•,•  «>  «,  i  >  -  t  i  >>  •  ,1 

rMO 


3 


VITA 


Dennis  James  Rensel  was  born  on  15  May  1950  in  Cleveland,  Ohio.  He 
graduated  from  Padua  Franciscan  High  School  in  Parma,  Ohio  in  1968  and 
attended  Case  Institute  of  Technology  of  Case  Western  Reserve  University 
for  one  year  prior  to  entering  the  United  States  Air  Force  Academy  in 
June  1969.  He  majored  in  both  Electrical  Engineering  and  Computer  Science. 

He  received  his  commission  and  B.S.E.E.  on  6  June  1973.  Upon  graduation 
he  attended  the  Communication-Electronics  Course  at  Keesler  AFB,  Mississippi 
from  August  1973  to  March  1974.  He  was  next  assigned  to  the  4754th  Radar 
Evaluation  Squadron,  Hill  AFB,  Utah,  where  he  served  as  a  Radar  Evaluation 
Officer.  Early  in  1976,  Lt  Rensel  attended  the  Air  Force  Communication 
Service  Wideband  Communication  Link  Evaluations  Course.  His  next  assign¬ 
ment  was  as  a  detachment  commander  of  a  remote  communication  site  in 
southern  Turkey,  TUSLOG,  Detachment  187.  One  year  later  he  moved  to  the 
European  Headquarters  of  Defense  Communication  Agency  (DCA)  as  a  Communi¬ 
cation  Engineer  working  in  Quality  Assurance  and  the  Transmission  Engineering 
Branches.  Elis  responsibilities  included  the  managing  of  several  joint 
service  installations  and  upgrade  programs.  In  June  1979,  Captain  Rensel 
was  reassigned  to  the  School  of  Engineering,  Air  Force  Institute  of 
Technology,  where  he  is  now  pursuing  an  M.S.E.E.  with  the  emphasis  in 
Communications  and  Digital  Systems.  He  is  a  member  of  Eta  Kappa  Nu. 

Permanent  Address:  3200  Hilltop  Drive 

Parma,  Ohio  44134 


132 


SECURITY  CLASSIFICATION  of  This  PAGE  'When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


1 .  REPORT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3  RECIPIENT’S  CATALOG  NUMBER 

Aw  x 


AFIT/GE/EE/80D-37  * 


4.  TITLE  (and  Subtitle) 

ERROR-CORRECTING  CODE  PERFORMANCE  IN  A  MODELED 
ECM  ENVIRONMENT 


au  Tho^j.i 

Dennis  J.  Rensel 


PERFORMING  organ  ration  NAME  ANG  ADDRESS 

Air  Force  Institute  of  Technology  (AFIT-EN) 
Wright-Patterson  AFB  OH  45433 


It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


776 


5.  TYPE  OF  REPORT  &  PERIOD  rovERED 


MS  Thesis 


6.  PERFORMING  OIG  REPORT  N.V9I.1 


8.  CONTRACT  or  gran  T  NUM0F  R 


io.  program  element  f»o;e:* 

AREA  &  WORK  UNIT  N  J  M  5  E  Ri 


12.  REPORT  DATE 

December  1980 


13  NUMBER  OF  PAGES 


14.  MONITORING  AGENCY  NAME  A  ADDRESSflf  different  from  Controlling  Office)  IS  SECURITY  CLASS  (of  th it  report. 

UNCLASSIFIED 

IS*  DECLASSIFICATION  OOWNGRaD'N.: 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (o I  this  Report! 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20.  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


Approved  for  public  release;] IAW  AFR  190-17 

4u>Aa^cC. 

FREDRIC  C.  LYNCH ,QMaj or. 


FREDRIC  C.  LYNCH^Ma  j  or ,  USAF  16  JUN  1391 

Director  of  Public  Affairs 


19.  KEY  WOROS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


Additive  White  Gaussian  Noise 
Probability  of  Error 
Viterbi  Decoding  Algorithm 
Dual-Three  Encoder 


8-ary  FSK 
Bit  Error  Rate 


^6.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  identify  by  block  number) 

Electronic  countermeasures  (ECM),  such  as  broadband,  switched  broadband,  contin¬ 
uous  wave  and  multitone  jamming,  are  designed  to  disrupt  communication  systems. 
The  subject  of  this  thesis  is  the  effect  of  these  ECM  on  the  performance  of  error 
correcting  codes  in  an  additive  Gaussian  noise  channel.  The  channel  is  modeled 
with  two  convolutional  coders  having  a  rate  of  h  and  constraint  lengths  of  2  and 
6,  an  8-ary  FSK  modulator  and  demodulator,  and  a  decoder  based  unon  the  V 1 1 o r u i 
algorithm.  Two  measures  of  effectiveness  are  the  probability  of  hit  errors  dm. 
for  the  waveform  channel  and  the  bit  error  rate  (BER)  for  the  overall  channel 


dd 


1473  EDITION  OF  1  NOV  65  IS  OBSOLETE 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  'M'S  PAGE  (When  On-'i 


_ UNCLASS  IK  l E1) _ 

ICCUWITY  CLASSIFICATION  OF  THIS  PAOEflWnn  Dmtt  Cnfrrdl 

performance.  It  is  shown  that  both  and  BER  are  increased  by  the  previouslyi 

mentioned  ECM.  To  aid  in  evaluating  the  performance  of  these  coders,  a  simula¬ 
tion  program  was  written.  In  the  program,  the  transmission  signal,  the  channel 
noise,  and  the  lamming  signal  are  independent  of  each  other.  The  basic  conclu¬ 
sions  for  these  two  convolutional  codes  are  that  the  dual-three  encoder  with  the 
longer  codeword  performs  better  in  the  no  jamming  or  in  a  Gaussian  jamming 
environment,  while  the  CC1  encoder  with  the  shorter  codeword  performs  better 
in  a  CW  jamming  environment.  ^ 

NOTE:  This  abstract  is  limited  to  200  words  maximum.  | 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  »C.F' *n- 


A-  -IT 


